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Abstract. Systems described by an 0{n) symmetrical 4> A Hamiltonian are considered 
in a d-dimcnsional film geometry at their bulk critical points. A detailed 
renormalization-group (RG) study of the critical Casimir forces induced between the 
film's boundary planes Q3j, j = 1,2, by thermal fluctuations is presented for the 
case where the 0(n) symmetry remains unbroken by the surfaces. The boundary 
planes are assumed to cause short-ranged disturbances of the interactions that can 
be modelled by standard surface contributions oc <fi 2 corresponding to subcritical 
or critical enhancement of the surface interactions. This translates into mesoscopic 
boundary conditions of the generic symmetry-preserving Robin type d n <p = Cj<p. RG- 
improved perturbation theory and Abel-Plana techniques are used to compute the 
L-dependent part / res of the reduced excess free energy per film area A — > oo to two- 
loop order. When d < 4, it takes the scaling form / rcs w D{c\L^/ v , c 2 £*/ !/ )/£ d ~ 1 as 
L — > oo, where Cj are scaling fields associated with the surface-enhancement variables 
Q, while $ is a standard surface crossover exponent. The scaling function £)(ci,c 2 ) 
and its analogue 2?(ci,c 2 ) for the Casimir force T = —df Yes /dL are determined via 
expansion in e = 4 — d and extrapolated to d = 3 dimensions. In the special case 
ci = c 2 = 0, the expansion becomes fractional. Consistency with the known fractional 
expansions of D(0, 0) and V(0, 0) to order e 3 / 2 is achieved by appropriate reorganisation 
of RG- improved perturbation theory. For appropriate choices of c\ and c 2 , the Casimir 
forces can have either sign. Furthermore, crossovers from attraction to repulsion and 
vice versa may occur as L increases. 
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1. Introduction 

When media exhibiting low-energy fluctuations are confined by walls or interfaces, or 
macroscopic objects are immersed in them, their fluctuation spectrum changes [1, 2, 3]. 
This may induce long-range effective forces between such objects and boundaries. 
A much studied example of such fluctuation-induced forces are the Casimir forces 
between a pair of parallel and grounded metallic plates produced by their influence 
on the quantum vacuum fluctuations of the electromagnetic field [4]. These quantum 
electrodynamics (QED) Casimir forces are known to be very weak unless the separation 
of the bodies between which they act becomes very small. The latter condition is met in 
micro-electromechanical devices (MEMS). As has recently become clear, Casimir forces 
of this kind must be taken into account in the design of such systems. Since they tend to 
be attractive for simple geometries, they may impair the functioning of MEMS, causing 
stiction [5, 6, 7, 8]. Owing to the interest in — and technological importance of — such 
small-scale systems, this insight has triggered considerable new activity in the study of 
QED Casimir forces. Crucial issues are the knowledge and control of their strength, 
sign, and geometry dependence [9]. 
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Another example of fluctuation-induced forces that has attracted a great deal of 
attention recently are the so-called thermodynamic Casimir force caused by thermal 
fluctuations near critical points [2, 10, 11]. Predicted decades ago [12], they were verified 
experimentally first in an indirect way through the thinning of 4 He wetting layers near 
the lambda transition [13] and their effects on wetting films of binary liquid mixtures 
[14, 15]. Subsequently they were measured directly in binary fluid mixtures near their 
critical points of demixing [16, 17]. 

In this paper we shall be concerned with thermodynamic Casimir forces in systems 
with a film geometry bounded by two planar free surfaces at a finite distance L. We 
will consider systems that can be modelled by an n-component |0| 4 Hamiltonian on 
sufficiently long length scales and study the Casimir forces directly at the bulk critical 
point. The internal 0(n) symmetry (or Z 2 symmetry if n — 1) is presumed to be neither 
explicitly broken by surface contributions to the Hamiltonian nor spontaneously for all 
temperatures T > T c oo , the bulk critical temperature. In other words, no ordered phase 
can occur for finite L when T > T Cj00 . Choosing a generic kind of such non-symmetry- 
breaking conditions, we will investigate the critical Casimir forces and demonstrate two 
important features of them. First, they can have either sign, depending on properties 
of the two surfaces — an expected result since attractive and repulsive critical Casimir 
forces were found a long time ago for certain boundary conditions that ought to apply 
on long length scales [18, 19, 20]. Second, for appropriate surface properties, crossovers 
from attraction to repulsion of the Casimir forces and vice versa can occur as a function 
of film thickness L. A brief account of parts of our work was given in [21]. The purpose 
of the present paper is to provide details of the calculations and further results. 

One motivation for our work is the obvious potential for cross-fertilisation between 
the fields of thermodynamic and QED Casimir effects. Since both types of effects share 
a number of characteristic features, mutual benefits may be expected. On the other 
hand, one must be aware of certain essential differences. Common to both types is 
that they exhibit universal properties and depend on gross properties of the fluctuating 
media and confining objects, as well as on their shapes and geometry. Two important 
differences are: 

(i) When studying QED Casimir forces, one frequently gets away with the analysis 
of effective Gaussian theories in which the coupling between the electromagnetic 
field and matter is accounted for by proper choices of boundary conditions. This 
usually holds even in the case of polarisable and magnetisable media where material 
properties also enter via dielectric and permeability functions. By contrast, 
adequate treatment of critical Casimir forces normally involve the study of non- 
Gaussian theories such as 4 theories or corresponding lattice models (e.g., Ising 
models) in bounded geometries. 

(ii) In studies of QED Casimir forces such as Casimir's original work [4], the 
electromagnetic fields are taken to have zero averages. Hence the Casimir forces 
are entirely due to fluctuations — if there were no fluctuations, there would be 
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no Casimir force. However, in the case of thermodynamic Casimir forces, the 
order-parameter densities 4>(x) have nonzero averages in ordered phases, where the 
order may either result from spontaneous symmetry breaking or else be imposed 
by symmetry-breaking bulk or boundary fields. If the medium undergoes a phase 
transition from a disordered high-temperature to an ordered low-temperature phase, 
then the order-parameter profile (<j>(x)} does not vanish in the latter. Such a 
nontrivial profile contributes to the size-dependent part of the free energy and 
hence, to the Casimir force. Consequently, a Casimir force will be obtained 
already at the level of Landau theory — i.e., in the absence of fluctuations - 
whenever it yields non-vanishing order-parameter profiles. Beyond Landau theory, 
the thermodynamic Casimir force will therefore consist of a part coming from a 
non-fluctuating background and a superimposed fluctuation-induced contribution. 

The last statement applies in particular to confined binary fluid mixtures [17, 
22, 23]. These are known to generically involve symmetry-breaking boundary fields. 
Therefore, the thermodynamic Casimir forces they yield will normally have contributions 
from non-fluctuating backgrounds on either side of the order- disorder transitions. Since 
even at the bulk critical temperature, {(f>(x)) is not expected to vanish, this holds there 
as well. 

It is well known that the boundary fields needed to describe binary fluid mixtures in 
contact with walls can have either sign, depending on which one of the two components 
gets preferentially adsorbed locally at the wall [10, 12, 22, 23, 24]. Furthermore, there 
is clear evidence based on theoretical and experimental work as well as on Monte Carlo 
simulations [10, 14, 15, 16, 25, 26, 27, 28, 29, 30] that the critical Casimir force in 
binary fluid mixtures confined between two planar walls can be attractive or repulsive. 
Corresponding crossovers must occur and were investigated at the level of mean-field 
theory [31]. 

Casimir forces in binary fluid mixtures are particularly well suited for direct 
experimental measurements. However, from our above remarks it is clear that their 
analogy with QED Casimir forces is limited: owing to the generic presence of a non- 
vanishing order parameter profile at, above and below the critical temperature, only 
a part of them is fluctuation induced. Since we focus here on the case in which the 
symmetry <fi — > —0 remains unbroken at T coo , we are dealing with critical Casimir 
forces that are entirely fluctuation induced and hence are more akin to QED Casimir 
forces, albeit due to thermal, rather than quantum, fluctuations. 

The rest of this paper is organised as follows. In section 2.1 the model is specified 
and the boundary conditions are recapitulated. Next, the free propagator is determined 
in section 2.2 for general values of the surface interaction constants, its eigenfunction 
representation is given and the transcendental equations derived whose solutions yield 
its eigenvalues. In section 2.3, many-point cumulant functions are introduced, their 
renormalization recapitulated and their RG equations presented. These considerations 
are extended to the free energy and Casimir force in section 3. In the remainder of 
this section details of our approach are explained and our analytical results for the 
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scaling functions of the L-dependent part of the free-energy density and the Casimir 
force are given. In section 3.5 it is shown how RG-improved perturbation theory can 
be reorganised to achieve consistency with the fractional series expansions in e one 
encounters in the special case of two critically enhanced surface planes. Section 4 
discusses our results, deals with the issue of how to extrapolate them to d = 3 dimensions 
and presents such extrapolation results for the scaling functions of the L-dependent 
part of the excess free-energy density and the Casimir force. Section 5 contains a brief 
summary and concluding remarks. Finally, there are 5 appendices describing various 
calculational details. 

2. Background 

2.1. Model and boundary conditions 

We consider a continuum model for a real-valued n-component order-parameter field 
<f){x) = (<f) a (x),a = l,...,n) defined on the <i-dimensional slab 2J = IR d_1 x [0, L\. 
We write the position vector as x = (xi, . . . , xj) = (y,z), where y = (yi, . . . , ya-i) 
denotes the d—1 Cartesian coordinates yj = Xj along the slab and z = Xd the remaining 
one across it. The thickness of the slab, L, is taken to be finite. We choose periodic 
boundary conditions along all y-directions so that the boundary <92J of 93 consists of a 
pair of (d— l)-dimensional planes 55 1 and 25 2 at z = and z = L. Assuming the absence 
of long-range interactions, we can choose the Hamiltonian to have the local form 

2 

-H = [d d x C h (x) + d'-'y Cj(y) , (2.1) 

where £5 and £,=1,2 are local bulk and surface densities depending on (f>(x) and its 
derivatives. 

In accordance with our considerations in the introduction we choose these densities 
to be 0{n) symmetric. For the bulk density we make the standard choice 

S 7- 7/ 

C h (x) = -[Vcj>{x)f + -4?{x) + ^|0 2 (*)| 2 (2.2) 

where (V0) 2 is the usual short-hand for J^ai^^a) 2 - ^ nas recently been emphasised 
that lattice systems lacking cubic symmetry, such as systems involving monoclinic or 
triclinic lattices, lead to 4 bulk densities with a gradient-square term of the generalised 
form ^g l i(di(f)) ■ dj<j> where <9j stands for the derivative d/dx l and is a non-diagonal, 
position-independent matrix [32]. We refrain from working with such generalised models 
because the matrix g l i is nothing but a constant metric [33]. The geometric effects it has 
can be absorbed by a proper choice of variables which of course enters the way lattice 
quantities are related to those of the standard bulk 4 model with action density (2.2). 
For details the reader might want to consult reference [33, p. 14-17]. 
The surface densities are given by 

£j(r) = f 2 M. (2.3) 
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The interaction constants c\ and 62 of the two boundary planes are allowed to differ. 
They are known to reflect the weakening or enhancement of the local pair interactions 
at 53 1 and 532, respectively [22, 23]. In semi-infinite systems bounded by a plane 53i, a 
threshold value c sp of the enhancement variable t\ exists such that a phase with long- 
range surface order appears (in sufficiently high dimensions) in a temperature regime 
above T Cj00 when the enhancement — 5c\ = c sp — ci is positive (is "supercritical"). The 
transitions that occur at T c oo in such semi-infinite systems are called ordinary, special 
and extraordinary, depending on whether the enhancement 5c\ is subcritical {8c\ > 0), 
critical (Sci = 0) or supercritical (Sc 1 < 0). Analogous statements hold for semi- infinite 
systems bounded by the plane 532 with surface enhancement 62- 

In order to rule out spontaneous symmetry breakdown at T = T Cj00 for large L, we 
require that both enhancement variables 5cj = Cj — c sp satisfy the condition Scj > 0, i.e. 
are non-supercritical. 

From the boundary contribution to the classical equation of motion we get the 
boundary conditions of Landau theory 

d n <j)(x) = Cj4>(x) for x e 53 i? (2.4) 

where d n is the inner normal derivative on <95J = 53 1 U532- Beyond Landau theory these 
boundary conditions still hold in an operator sense (inside of averages) [22, 23, 34, 35, 36]. 



2.2. Free propagator 

A quantity of central importance for the renormalization-group (RG) improved 
perturbation theory we are going to use below is the free propagator Gl at r = 0. It is 
given by the operator inverse Gl = (— A) -1 subject to the boundary conditions (2.4), 
where A = V 2 is the Laplacian. Let \m) be the complete set of eigenfunctions of the 
operator — d 2 z = —d 2 /dz 2 on the interval [0, L\ satisfying the boundary conditions 

(d z -c 1 )\m)\ z=0 =0, (2.5) 
(d z + c 2 )\m)\ z=L = 0. (2.6) 

That is, we have 

-d 2 z \m) = k 2 m \m), (2.7) 
and the eigenstates fulfil the orthogonality and completeness relations 

(m\m) = 8 mm > (2.8) 

and 

00 

](z\m)(m\z') = 6{z-z'), (2.9) 



m=l 



where the label m = 1,2, . . . , 00 enumerates the eigenvalues by size, starting with 
the smallest one. 
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To exploit the translation invariance of the system along the direction parallel 
to the boundary planes 93 j, we Fourier transform with respect to the difference of y 
coordinates, writing 

G L (x,x')= ^ G L { P] z,z')e ip <y'-y\ (2.10) 



(2.11) 



where we have introduced the notation 

"(d-l) r d d - 1 p 

= L-! (27T) d - 1- 

We thus arrive at the spectral representation 

a. , i \ (z\m)(m\z') 

G'to'^-^-ttk 1 - (2 - 12) 

m=l 

Since our assumption that both enhancement variables 6j are non-supercritical 
implies that the associated renormalized variables Cj (whose definition will be recalled 
in the following section 2.3) are nonnegative, we shall need the eigenvalues k^ and 
eigenfunctions \m) only for Cj > 0. In this case, there are infinitely many eigenvalues 
for finite L, all of which are nonnegative and non-degenerate (see reference [37, 
Appendix A] and below). A zero eigenvalue k\ — occurs for arbitrary L e (0, oo) only 
when Ci = c 2 = 0. 

Let us briefly summarise the relevant properties of the eigensystem we 
shall need in our subsequent analysis. On dimensional grounds the eigenfunctions must 
have the form 

(z\m) = L~ 1/2 v m (zjL\hL, c 2 L). (2.13) 
It is convenient to introduce the dimensionless variables 

( = z/L, n m = Lk m , Cj=CjL. (2.14) 
Setting L — 1 in equation (2.13) one concludes that the functions v m can be written as 

Vm((\Cl,C 2 ) = \]l 2 COs{K m C-$ m ). (2.15) 

The analogue of the boundary condition (2.5) 

^(0|Ci,C 2 )=Ci7; ro (0|Ci,C2) (2.16) 
implies that the phase shift ^ m is given by 

tani? m = Ci/K m = cxjk m . (2.17) 
We can choose it such that < -d m < ir/2, which leads to 

sin^ m =^= 1 cos $ m= = m / Cl . (2.18) 

A/l+(/t m /Cl) 2 v 1 + ( / W c i) 

The limiting values $ m = and <d m = n/2 are obtained when c± — > and c.\ — >■ oo, 
respectively. 
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Using these results along with equation (2.15), the normalisation constant A m can 
be computed in a straightforward fashion. One obtains 



A, 



1+ Cl ° 2 



(2.19) 



To determine the spectrum {n^}, we use the analogue of the boundary condition (2.6), 
v' m (l\C 1 ,C 2 ) = -C 2 v m (l\C 1 ,C 2 ). (2.20) 

This yields the transcendental equation 

R Cl ,c^ m ) = (2.21) 

for K m , with 

R Ci ,c 2 (k) = (C 1 +C 2 )kcos(k) + (dC 2 - K 2 )sinK. (2.22) 

Equations (2.5)-(2.7) specify a regular Sturm-Liouville problem for which the 
following mathematical properties are known (see e.g. references [38], [39, Kap. IX], 
[40, Kap. IV.3] and [41, chapters 8, 9]): (i) The eigenvalues K 1 m {C\,C 2 ) are real, non- 
degenerate, countable and accumulate only at oo. (ii) They can be ordered such 
that < k 1 ^ for m < m! . There is a smallest eigenvalue k\ but no largest one, 
i.e. K m — > oo as m — > oo. (iii) The eigenfunctions v m corresponding to different 
eigenvalues are orthogonal with respect to the standard L 2 scalar product in L 2 (0, 1).| 
(iv) The normalised eigenfunctions are a complete orthonormal set (basis) in the 
Hilbert space L 2 (0,1). (v) When C\ and C 2 > 0, we can use the Robin boundary 
conditions (2.5) and (2.6) in conjunction with the fact that the eigenfunctions vanish 
only for Dirichlet boundary conditions at the boundary planes Q3j to conclude that 
v(dv/d()\^zl = — C 2 v 2 \^ = i — CiV 2 \^ = o < 0. Fulfilment of this condition guarantees (by 
theorem 7 of reference [40, p. 234]) that k\ and hence all strictly positive. 

The non-degeneracy of the eigenvalues be verified explicitly by showing 

that the positive zeros K m of -Rci,c 2 ( K ) are simple. To this end, one can compute the 
derivative R'c lt c 2 ( K ) = <9k-Rci,c 2 ( k ) at k = k to and express the trigonometric functions 
cosn m and sin/t m in terms of K m , C\ and C 2 using equations (2.21) and (2.17) along 
with the analogue of the latter implied by the boundary condition at ( — z/L — 1. This 
gives 

«4 + + c 2 + c! + cp + c 1 c 2 (c 1 + c 2 + CiC 2 ) 
V« + + c 2 2 ) 

which is nonzero for all C\ > 0, C 2 > and n m > 0. 

For the special choices (Ci, C 2 ) = (0, 0), (0, oo), (oo, 0) and (oo, oo), equation (2.21) 
can be solved explicitly for the k^- The respective results for Neumann-Neumann 
(N, N) Neumann-Dirichlet (N, D), Dirichlet-Neumann (D,N), and Dirichlet-Dirichlet 
(D, D) boundary conditions can be looked up in Appendix A of reference [20]. However, 
for general values of the Cj, the zeros of the transcendental equation (2.21) cannot 

\ In the general case of Sturm-Liouville problems, the scalar product involves a weight function [39]. 
This is unity in our case. 



r' c c ( K ) = (-i) m m mK 1 T^gJLrLZ: ^ T t^t uiu 2 ; ^ 23 ^ 
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be determined in closed analytical form. This is a familiar difficulty, which makes the 
evaluation of the single and double mode sums J2 m one encounters in the calculation 
of one and two-loop Feynman graphs of the free energy a nontrivial problem. We shall 
turn to this issue in section 3. 

Note also that in some of the above formulae we have tacitly assumed that C\ and 
C 2 are both positive. However, these equations remain valid, firstly, for all m when 
either one of the Cj vanishes while the other remains positive, and secondly, for all 
m > 1 when C\ = C 2 = 0, because the respective n m then all approach positive values 
in the corresponding limits. Hence one can simply set C\ or C 2 (or both) to zero in the 
respective equations. The case of m — 1 with Ci,C 2 — > is special in that K\ — > 0. It 
is an easy matter to check that the correct limiting eigenfunction -u(£|0,0) = 1 results, 
independent of the order in which the limits C\ — > and C 2 — > are taken. 

It is also not difficult to determine how the eigenvalues k^(Ci,C 2 ) approach the 
(N, N)-values K m (0, 0) — (m — l) 2 vr 2 , m — 1, 2, . . . , 00, as both Cj approach zero. One 
can either determine the terms linear in Cj from equation (2.21) or else perform Rayleigh- 
Schrodinger perturbation theory in the potential C\ 5(() + C 2 5(1 — C) to see that 

«4(d,C 2 ) « (m-l)V + (2-5 m , 1 )(C 1 + C 2 ), m = l,2,...,oo.(2.24) 

Ci,C2— >o 

2.3. Many-point cumulants, their renormalization and RG equations 

We proceed by providing some of the necessary background on the renormalization of the 
theory defined above. § Let us begin by recalling the bulk and boundary counter-terms 
required to absorb the ultraviolet (UV) singularities of many-point cumulant functions. 
To this end we consider the (A^+M 1 +M 2 )-point cumulant functions 

r (N;M u M 2 ) , \ 
G a 1 ,...,/3M 1+ M 2 ^ 1 ' • • • >2/M 1+ M 2 J 

IN Mi M 1 +M 2 \ cum (2.25) 

\j=l k=l «=Mi+l / 

involving N fields 4> aj (xj) located at points in the interior of 2J and Mi+M 2 fields 4>p k 
on the boundary planes Q3i and Q3 2 . As indicated, we take the first Mi boundary points 
to lie on 03 1 and the remaining ones on 53 2 - 

Equations (2.12)-(2.17) imply that the p transform of the free propagator for 
general nonnegative values of b\ and c 2 can be written as 

G l (p;z,z'\t,ci,c 2 ) 

1 \- [km cos(k m z) + 6i sin(k m z)][k m cos(k m z') + ci sin(k m z')] (2.26) 
= (p 2 + r + kl)(cj + kl) 

Its explicit form is known [33, 42]. It is given by the expression into which the right-hand 
side of equation (B5) of reference [33] turns upon replacement of its by 

Xp = VV + r, (2.27) 

§ A more extensive discussion of this issue can be found in reference [22] . 



10 



i.e. 



G l (p;z,z'\t,c 1 ,c 2 ) 

q. , , [c\ smh(x p z) + x p cosh(x p 2;)] [c 2 sinh[x p (L — z')\ + x p cosh[x p (L — z')] 
Xp(6i + c 2 ) cosh(xpL) + Xp(x^ + cic 2 ) sinh(xpL) 
+ (2 <->■ z'). (2.28) 

Both expressions (2.26) and (2.28) are exceedingly difficult to work with. To 
understand the nature of the UV singularities of the theory, it is better to treat the 
boundary terms Cj as interactions and work with the free propagator for C\ = c 2 = 0. 
The latter is nothing but the free propagator G^ N satisfying Neumann (N) boundary 
conditions on both boundary planes. It can be written as [22, section IV. A] 

Gr(p;z,2>) = G L (p;z,2>\T,0,0) 

(2.29) 



= [<?b(p, z~z' + J2L) + G h (p, z + z' + j2L)], 



where 



is the free bulk propagator for f = in the pz-representation. As expounded in 
reference [22], the j = contribution from the first term in square brackets in 
equation (2.29) is the origin of the familiar primitive bulk UV singularities. The j = 
and j = — 1 contributions from the second term in square brackets are singular at 
coinciding points on 23 1 and 93 2 , respectively. The additional UV singularities they 
produce can be absorbed by counter-terms with support on these boundary planes, 
i.e. the surface counter-terms required for the corresponding semi-infinite systems 
[22, 35, 36]. All other contributions in the square brackets do not diverge at coinciding 
points and hence do not give rise to additional UV singularities. The upshot is that the 
usual bulk reparametrizations 

4> = zJ /2 0r, 

t - t c ,oo = St = Z T fi 2 r, (2.31) 
uN d = fi 4 ~ d Z u u, 

of the field, temperature variable f and interaction constant u in conjunction with the 
reparametrizations 

Cj ^sp = SCj f^Z c Cj, 

of the surface enhancement variables 6j and the boundary operators (fy Bl {y) = <f>(z,0) 
and <B2 (?/) = 4>(z, L) suffice to absorb the UV singularities of the functions G^ N ' Ml ' M2 \ 
Following reference [43], we choose the factor that is absorbed in the renormalized 
interaction constant as 



N = 2r(3-d/2) 1 



(rf-2)(4vr) d / 2 16tt 2 



l+ 1 -^: ln(47r) e + Q(e 2 ) 



(2.33) 
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where e = 4 — d and 7^ is the Euler-Mascheroni constant. This choice ensures that 
the bulk renormalization factors Z^, Z T and Z u as well as the surface renormalization 
factors Z\ and Z C) when determined by minimal subtraction of poles in e, up to second 
order in u reduce to the two-loop results given in references [22], [34], [35], [36], [44] 
and used in Krech and Dietrich's work [19, 20]. The quantity f Cj00 is the critical bulk 
value of f . Further, c sp is the critical enhancement value associated with the special 
surface transition. In a theory whose UV singularities are regularised by means of a 
large-momentum cut-off A, these quantities diverge ~ A 2 and ~ A, respectively. In our 
calculations below we shall utilise dimensional regularisation. 

Note also that the surface renormalization factors Z\ and Z c depend exclusively on 
u and e but not on Cj (nor on L) when fixed by the requirement that the UV poles in e 
be minimally subtracted. In our calculations of the L-dependent part of the free energy 
at the bulk critical point and the critical Casimir force to be described in section 3 we 
shall need Z c merely to first order in u. We therefore quote the result 

Z c = l + — —u + 0{u 2 ) (2.34) 

for convenience. 

Upon introducing the renormalized cumulants 

^....aJ^r^i, r, Cl , c 2 , L, n) 

7 -N/2 (7 y ) -(M 1 +M 2 )/2 MN;Ml,M2) / .iii ZT\ ^ ' ' 

we can now derive RG equations for these functions by varying /i at fixed bare parameters 
u, f, ci, c 2 , and L. Let us define the exponent functions 

ri g = fid^Zg, g = u,T,c 1 ,c 2 ,l,<f>, ( 2 - 36 ) 
the beta function 

(3 u (u) = fid^\ Q u = -[e + r} u (u)]u, (2.37) 

(where d^\ indicates a /x-derivative at fixed bare parameters u, f, ci, c 2 , L) and the 
operator 

= ^ + f3 u du - (2 + r]r)rd T - (1 + r] c )( Cl d Cl + c 2 d C2 ). (2.38) 
Then the RG equations can be written as 



' N Mi + M 2 . 

T>„ + — rjcp + g ^ + W 



They are completely analogous to those for the renormalized cumulants of the respective 
semi-infinite geometries (L = 00 with either M 2 = or Mi = 0). We have given them 
here for general r, even though we will restrict ourselves to the critical case r = in 
our calculations in section 3. 
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3. Free energy and Casimir force 

3.1. Definitions 

We next turn to the free energy and the Casimir force. Introducing the partition function 
Z and the total free energy F via 



e -HM = e -F/k B T ^ (31) 



we define the reduced free-energy density per hyper-surface area A, 

F 

/(L;f,£,ci,c 2 ) = Jim tt^tt- (3-2) 

This quantity can be decomposed as 

f(L;r,u,c 1 ,c 2 ) = Lf h (r,u) + f s + f Tes (L; f, u, Ci, c 2 ), (3.3) 

where /b and / s are the bulk free-energy density per <i-dimensional volume LA — >■ oo 
and the excess free-energy density per hyper-surface area A — > oo, both of which are 
defined by appropriate thermodynamic limits (see e.g. references [22] and [45]). The 
former depends only on the bulk interaction constants, the latter additionally on the 
boundary interaction constants. It is a sum of the reduced surface excess free-energy 
densities f B j of the corresponding semi-infinite systems bounded by either < S 1 or Q3 2 : 

f s (T,U,C 1 ,C 2 ) = / S) i(f,«,Ci) + / s ,2(r,M,c 2 ). (3.4) 

The remaining term in the decomposition (3.3), the reduced residual free-energy density 
/ res , contains the full L-dependence of / — Lf^. The effective ("Casimir") force per 
hyper-surface to which it gives rise is 
HL;T,...) df Tes 

k B T =—dL> (3 ' 5) 
where the ellipsis stands for all other (bulk and surface) variables on which it depends. || 



3.2. Renormalization of the residual free energy and Casimir force 

Inspection of the perturbation series of the reduced free energy F//c B T reveals that the 
bulk reparametrizations (2.31) and surface reparametrizations (2.32) are not sufficient 
to absorb all UV singularities. This is because both the bulk free-energy density /b as 
well as the surface free-energy densities f s j have additional primitive UV singularities 
[22]. To cure these, additive bulk and surface counter-terms are required. Since these 
counter-terms must merely absorb primitive UV singularities of /b and f s j, respectively, 
they can be chosen independent of L. One convenient way to fix them is by subtracting 
from /b the Taylor expansion in 5f, and from f s j that in St and 56j about non- 
vanishing reference values to the appropriate orders (see e.g. [22, section III. C. 12] or 
[43, section II. E]). Owing to the L- independence of these additive counter-terms, their 

|| Since the right-hand side of equation (3.5) depends on the variables f, u, b\ and 62, so does its 
left-hand side. 
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contributions cancel out in the difference of free-energy densities that / res involves. Using 
the reparametrizations (2.31) and (2.32) therefore gives us the UV-finite renormalized 
residual free-energy density 

f reS;R (L;T,u,c 1 ,c 2 ,ii) = f Tes (L;T,u,c 1 ,c 2 ). (3.6) 
Hence this quantity satisfies a homogeneous RG equation, namely 

?Vres,R = °- ( 3 - 7 ) 

The latter can be solved in a standard fashion by means of characteristics. Let us 
define the running variables u(£) and Cj(£) as solutions to the initial value problems 

e±u(£) = p u [u(e)], u(l)=u, (3.8) 

l^fi® = "{I + VcWmW, cj(l) = c r (3.9) 

Let u* be the infrared-stable zero of the beta function f3 u for d = 4 — e < 4. We shall 
need its e expansion 

3e 9(3n + 14)e 2 t 
n + 8 + (n + 8) 
below only to 0(e). Solving the flow equations (3.9) for Cj yields 

(3.11) 



u* = Z^+ ,",oT +°( £3 ) (3- 10 ) 



E c (u,u) = exp j- J du'[l + r] c (u) -y c ]| 



where $ is the surface crossover exponent whose e expansion is known to 0(e 2 ) [35, 36]. 
Since we shall need the e expansion of the RG eigenexponent y c to first order in e, we 
recall the result 

l/c = l-^ e + 0(e 2 ). (3.12) 

In the large- length-scale limit £ — > 0, the running variables u and Cj behave as 
u{£) ~ -u* and 

c 3 (£)^ o E:(u)£-^c 3 , E* c (u)=E c (u*,u) : (3.13) 

respectively, where E*(u) is a non-universal amplitude. We now set r = 0, solve the 
RG equation (3.7), choose the flow parameter as l\xL = 1 and use the above limiting 
expressions for u and Cj. This gives the asymptotic behaviour 

f ies , R (L;0,u,c 1 ,c 2 ,Li)/n w ^^zr" (3.14) 



with 



where 



D(d, c 2 ) = / res (l; 0, M* , Ci, c 2 , l)/n (3.15) 
c, = E:( M ) Cj (/iLr. (3.16) 
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The analogous scaling form of the reduced critical Casimir force per hyper-surface 
area follows by differentiation with respect to L. It reads 

J r {L;T c ^ 00 ,u,c 1 ,c 2 ,fi) 



nL- d V(c u c 2 ) 



with 



P(ci, c 2 ) = [d-l + y c {cid Cl + c 2 d C2 )]D(c u c 2 ). 



(3.17) 



(3.18) 



To appreciate these results, recall that the large- L form of the residual free-energy 
density at the bulk critical point T = T c oo , for given asymptotic large-scale boundary 
conditions BC, is conventionally written as 



f BC (T-T \ ~ A BC T-( d - 1 ) 

L^oo 



(3.19) 



where A BC is an L-independent universal, but BC dependent, number, called "Casimir 
amplitude". According to our result (3.15), a scaling function D appears for general 
values of c\ and c 2 in place of the Casimir amplitude — that is, the Casimir amplitude 
becomes scale- dependent. 

On the other hand, various Casimir amplitudes A BC can be recovered from the 
knowledge of _D(ci,C2). Since we assumed both surface enhancement variables Cj to be 
non-supercritical, and also ruled out the breaking of the 0(n) symmetry via boundary 
terms, the case of Robin boundary conditions we consider includes the four cases 
BC = (0,0), (0,sp), (sp, O) and (sp,sp) of large-scale boundary conditions that are 
associated with the fixed-point values Cj = oo and Cj = of the respective ordinary (O) 
and special (sp) fixed points of the semi-infinite systems bounded by *Bj. Hence we have 
A (o,o) /n =jD ( 00)00 ) ) 

A (o,s P )/ n = D ^ ) = A (sp,o)/ n = £)(o, oo), (3.20) 
A (s P ,s P )/ n = ^(0,0). 

The first two of these amplitudes, A^ 0,0 ) and A^°' sp \ are known to have expansions 
in integer powers of e. The leading two terms of these series expansions were determined 
in reference [20]. The corresponding results can be written as 

A (o ^/n = a + a 1 (n)e + O(e 2 ), (3.21) 



and 



with 



and 



A (0,sp) /n = 



-7a 



+ 



7T 



n + 2 4 In 2 



1024 V n + 8 



45 



7ax(n) 



a 



7T 



1440 



e + 0(e 2 ) (3.22) 



(3.23) 



ai(n) 



7T 



2880 



1 - 7s - ln(4vr) + 



2('(4) 5n + 2 
C(4) + 2 n + 8 



(3.24) 
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By contrast, the Casimir amplitude A^ sp ' sp ^ does not have a power-series expansion in 
e. As shown in references [43] and [46], the presence of the k\ = mode leads to 
a breakdown of the e expansion of A (sp,sp * ) and produces additional half-integer powers 
e^ 2 with j — 3, 5, . . ., which are modulated by powers of In e when j > 5. The expansion 
is known to order e 3 / 2 ; it reads [43, 46] 

A( sp ' sp 7n = a + ai(n) e + a 3/2 {n) e 3 / 2 + o(e 3/2 ) (3.25) 

with 

The e-expansion results for the scaling function D(ci,C2) we are going to derive 
below must be consistent with the series expansions (3.21)-(3.25) and equation (3.20). 
This requirement provides nontrivial checks for these results. 



3.3. Perturbation theory 

We next turn to the loop expansion 

oo 

f{L- f, u, c 1; c 2 ) = f [l] «» Ci, c 2 ) (3.27) 

of the reduced free-energy density and its computation for f = to two-loop order. The 
zero-loop contribution vanishes in the disordered phase. The next two terms can be 
written as 

/W (L; f, u, c!, c 2 ) = ^ Q = ^ Tr In (f - A) 

oo (3.28) 
= ££/ Mr + ^ + A:^ 

Z m=l 

and 

-•[21/ r o „ „ „ N — 1 C J ™(™ + 2)« /" d d X . . „ 2 

/ l2J (L;r,M, Cl ,c 2 ) = — W= - y — [Gifoslci.ca)] 



un(n + 2) ^ Z^" 1 ) /•( d - 1 ) A mi , mi , m2jm2 (c 1 L,c 2 L) 



4! mi £=X ^ (p? + + (p! + *2» + *) 

for general f > 0. Here we have introduced the functions 

rl 4 

"mi ,m2 ,ni3 ,m4 

(Ci,C 2 )= / dCMu roj .(C|Ci,C 2 ). (3.30) 

For the parameter values for which they are needed in equation (3.29), a lengthy but 
straightforward calculation yields 

^mi,mi,m2,m2 (fill ^2) 

1 V^mi ^rn 2 K m 2 p (<r,p) 1 x r (3.31) 

" 8 ^ (C? + <) (C 2 + <) (C 2 + «*J (C 2 + ^ + 4 Amidmi ' m2 



(3.29) 
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with the polynomials 

P^l = 2C 3 1 Cl(C 1 + C 2 + C 1 C 2 ), 

= 2[Cf + Cl + 2C 2 C 2 + 2dC 2 2 + {C\ + C 2 2 ) 2 ], 

^2 = 2, 

p£;g = 2C 1 c 2 (c 1 2 + cf)(d + c 2 + dC 2 ) = p#g, 

P^ 2 = 2{C l+ C 2 + Cl + Cl) = P^ 

P^g = 2C 1 C 2 (C 1 + C 2 + dC 2 ) = P c ( °g. (3.32) 

The p-integrals in the above equations (3.28) and (3.29) can be handled in a 
standard way using dimensional regularisation. The main difficulty one then is faced 
with is the calculation of the resulting single and double sums over the not explicitly 
known eigenvalues k^. We have done this by means of complex integration, modifying 
the Abel-Plana techniques described in reference [37] for our purposes.^ The technical 
details are given in Appendix A-Appendix B. Here we just present our results. 

As before, we use the notation to specify the /-loop term of a quantity Q. 
We choose the additive constant in the free-energy density such that the dimensionally 
regularised bulk free-energy density /b vanishes for f = 0. Hence 

i 1] (r = 0) = fl 2] (r = 0,u) = 0. (3.33) 

Our results for the surface free-energy densities read 



and 



^ =o ^) = - 2 (d -i)S: w ^ A (3 - 34) 

fM(T -One-)- n(n + 2) * r2 [ (3 ~^ )/2] g (3 35) 
} sj{T -U,u, Cj )- 2 4 , cos(d7r) _ 1 , 

where we have introduced the familiar quantity 

r(d) 

K d = 5(q-l) = 2 1 - d n- d ^ 2 /T(d/2). (3.36) 

Jq 

In order to write the one- and two-loop residual free-energy densities in a compact 
fashion, it is helpful to define the functions 



i _ (t_d)(t-c 2 ) e _ 2t 



(t + d)(t + C 



2) 



(3.37) 



Y <to) 1 (d + 2v \ C? +d ~* - C 2 - +d - 4 



2 /"°° f-l)" t2ff+(i-3 

YS> = --co S (*/2)/ o d' (f + Cl)2(i + c 2) 4 2 .-y- C? )( f2 - C |) < 3 ' 39 ' 

As we show in Appendix A. 2, it can also be derived from the result (2.26) for the free propagator. 
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and 



Ad) 

J C\,C2 



sin(<i7r) 



dtt 



2d- 



l _ (t-d)(t-c 2 ) e _ 2t 



{t + C 1 ){t + C 2 ) 



In terms of these, our results can be written as 



Ad) 
Lc\ ,Lc,2 



and 

f®(L-0 u & &,) - n(n + 2) fi r2 [( 3 - rf )/ 2 ] L 5- 2 , 

J ICS {L, U, U, Ci, C 2 J - - - ^y-i L 



/ j Lc\,Lc2 \ Lc\,Lc2 Lc\,Lc2 Lc\,Lci Lc\,Lc2 

a,p=0 



(3.40) 



(3.41) 



(3.42) 



3.4- Renormalized residual free- energy density and scaling functions 

The functions Xj££\ L and xj£'£ \ L have simple poles at d — 4, caused by the behaviour 
of the respective pre-factor of the integral in equation (3.38). These UV singularities 
imply that the two-loop term (3.42) is not regular at e = 0. It has a simple pole 
originating from the terms proportional to X^^\ 2L and X^\ 2L . One finds 

f®(L; 0, u, ci, c 2 ) n(n + 2)uL~'N d 



(c^ + c^^L 2 ^, 



) 



£-(d-i) Xi - 7 V ^ " C1 ^ 2 ^ (3.43) 

+nKU) + o( e °), 

where uiV rf = //w + 0(w 2 ). It is an easy matter to check that this pole gets cancelled 
by the contribution oc u/e one obtains from fies(L; 0, Ci, c 2 ) upon expressing the bare 
variables Cj in terms of their renormalized analogues Cj via equations (2.32) and (2.34). 
Substitution of the one- and two-loop terms (3.41) and (3.42) into equation (3.6) 
therefore yields indeed a UV- finite 0(u) result for the renormalized residual free-energy 
density f Tes , R {L; 0, u, c 1 ,c 2 , //), namely 



f TCS:R (L;0,u,c 1 ,c 2 ,n) 
nL-^- 1 ) 



--D (ciLfi, c 2 Lfi) + e 



1 - ^ E } D (ciL/jL, C 2 L/jL) 



1 [°° 

~ ItT 2 J dth ciL^,c 2 L^(t)t 2 ^t 



+ U 



n + 2 



1 

5^ [ ~~ 2 + ln C A 3 A)(ciL/i, c 2 L[i) 



(3.44) 



+ 



n + 2 
12vr 2 



EpM /(") i + OlVl 

ciLfi,C2Lfi ciL/i,C2Lfi ciLfi,C2Lfi I 1 W V U ' /! 



<r,p=0 



where 



D Q {C U C 2 ) = /Pl(l;0,Ci,C 2 )/n = ^ jf dtt 2 /ic lA (t) , 



(3.45) 
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and 



jto =--Y (4 ' a) 



dt- 



(-l)'f 



(T J-l + 2cr 



( f + Cl ) 2 (t + C 2 ) 2 e 2 * - (t 2 - CIW - C 2 2 ) ' 
To obtain the e expansion of the scaling function D, we set u = u* . Using 

n + 2 



(jjlL)* = imL 



1 



n + 8 



e ln(/iL) 



+ 0(e^ 



(3.46) 



(3.47) 



one sees that the result is consistent with the predicted scaling form (3.14) and yields 
the expansion 

D(d, c 2 ) = D (ci, c 2 ) + e£>i(ci, c 2 ) + 0(e 2 ) 



(3.48) 



with 

A(ci,c 2 ) 



1 - 



7£ - In 7r 



1 Z" 00 

j£> ( Cl ,c 2 )-— ^ dt/i CliC2 (t)t 2 lnt 



+ 



n + 2 
n + 8 



E 

j'=i 



1 

2 



ln(2 Cj 



Cj9 c .Do(ci,c 2 ) 



(3.49) 



+ 



i 2 

_2 / / Ci,C 2 tV Ci,C 2 J Cl,C 2 



4tt 2 



cr,p=0 

For the special values c, = and oo, the functions D and D\ can be analytically 
computed in a straightforward manner. One finds 

A)(oo, oo) = A)(0, 0) = — D (0, oo) = 



71 



Di(cx), oo) = £>i(0, 0) = ai(n), 



£>i(oo,0) 



7T 



1024 



n + 2 4 In 2 



1440' 
7ai(n) 



(3.50) 



s n + 8 45 y 8 ' 

where ai(n) is the coefficient defined in equation (3.24). These results confirm the 
validity of the relations (3.20) to first order in e. However, the present form (3.48) 
of our result does not yield the contribution ~ e 3 / 2 to D(0, 0) = A( ' sp )/^- The 
reason is that we assumed that c\ + c 2 > so that the lowest eigenvalue k\ of —d\ 
[cf. equation (2.7)] is strictly positive. Since k\ = when c\ = c 2 , there is a zero 
mode in the free propagator. As already mentioned, this causes a breakdown of the 
conventional RG-improved perturbation theory at T Cj00 [43, 46, 47]. In the following 
we will improve our results in such a manner that they fully comply with all of the 
small-e expansions of Casimir amplitudes given in equations (3.21)-(3.25), including 
the fractional expansion (3.25) of A( sp ' sp ) to 0(e 3 / 2 ). 



3.5. Modified RG-improved perturbation theory 

3.5.1. Formulation and results It is known from references [43], [46] and [47] how 
one can cope with the mentioned zero-mode problem one encounters at r = when 
c i = c 2 = 0: one can reorganise RG-improved perturbation theory such that it becomes 
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well-defined at r = 0. This suggests an obvious strategy to ensure consistency with the 
small-e expansions of A( sp - sp ) to e 3 / 2 . We should reorganise RG-improved perturbation 
theory for general c\ > and c 2 > in a way that reduces to the one used in previous 
work [43, 46] for the case c\ = c 2 = 0. 

This can be done as follows. We decompose the order-parameter field into an 
(m = l)-component and an orthogonal remainder, writing 

0(1/, z) = L- l / 2 ip(y) Vx(z/L) + if)(y, z), 

(3.51) 



m=2 



with 



• (to) 



& m) (v)= / d(v m (0<f>(yXL) (3.52) 



n 



and <p(y) = 4>^\y). The main difference compared to the previously studied cases of 
periodic and (sp, sp) BC [43, 46], as well as generalisations of the latter [47], is that 
the (m = l)-mode is not z-independent unless c\ = c 2 = 0. Thus the orthogonality 
of the ipa and the eigenfunction v± does no longer translate into the vanishing of the 
integral dz if>(y, z). The 4 term of the Hamiltonian therefore now generates also a 
vertex oc (if) • tp)(f 2 (which vanishes when c\ — C2 = 0), in addition to four-point vertices 
involving I = 0, 1, 3 and 4 fields ip a and 4 — / powers of ip a . 

Let us introduce fy, the free-energy density associated with the if)- field, by 

^ = -lnTr^e- ww . (3.53) 

Its loop expansion 

f 1 p(L;f J u,c u c 2 ) = /j 1] (L;f,M,Ci,c 2 ) + f [ ^ ] (L; r , u, c u c 2 ) + 0(3-loop) 

-1 1 (3.54) 

= ^"'^ >< + 0(3-loop) 



is analogous to that of /(L; f, u, ex, c%) specified in equation (3.27). The one- and two- 
loop terms and are given by the series in the second lines of equations (3.28) 
and (3.29) excluding the summands with m = 1 and m 1 ,m 2 = 1, respectively. In the 
corresponding Feynman graphs (second line), we depict free ^-propagators as dashed 
blue lines. 

We now set f = and define an effective (d — l)-dimensional field theory with 
Hamiltonian H e s[(p] by integrating out if): 

e -WH = ^ e -«M, (3.55) 

To determine "Hefrf^], we use perturbation theory. The contribution to the two-point 
vertex caused by the coupling between ip and if), to first order in u, originates from the 
graph 



1 

V / 



Sr L (ci,c 2 ) / d d ~W, (3-56) 
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where 

iTi(ci ' C2) — (3 - 57) 

It represents a local interaction corresponding to a shift 5tl of f. The Hamiltonian 
Hesi^p] becomes 

n cS [cp] = [ d^y \\{V v vY + \ [5t l {c x ,c 2 ) + kf\ ^ 

3 . 1 , (3.58) 

u 1 

+ ^ z A 1 , 1 , 1 , 1 (c 1 L,c 2 L)|^| 4 + 0(w 2 )j. 

It should be clear that beyond first order in u also nonlocal two-point and 2s-point 
vertices with s > 2 appear in T-L e s[ l p]- Two examples of such graphs are shown in 
figure 1. 



¥> <f 

i ii i 

+ % * y cp' 1 



x: :x 

, — s 

Figure 1. Examples of graphs producing nonlocal two- and four-point vertices of H c s . 



To obtain the free-energy density f(L; 0, u, &x, c 2 ) in the present modified 
perturbation scheme, we must add to 0, u, ci, c 2 ) the contribution associated with 
the y-field. We denote it as and define it via 

A = - In Tr v e- WcffM (3.59) 

by analogy with equation (3.53). Hence we have 

f(L; 0, u, ci, c 2 ) = U{L; 0, u, ci, c 2 ) + f$(L\ 0, u, Ci, c 2 ), (3.60) 

where /^(L; f = 0, u, Ci, c 2 ) is given by equation (3.54) and /^(L; 0, m, Ci, c 2 ) has the loop 
expansion 



f ip (L;0,u,c 1 ,c 2 




(3.61) 



= /W (L; 0, u, c 1 , c 2 ) + ff (L; 0, u, c ls c 2 ) + . . . . 

Here and /| ] are the one- and two- loop terms of a (d — l)-dimensional </> 4 theory 
with quadratic and y? 4 interaction constants SrL + kf and uAi^i/L, respectively. Their 
explicit expressions may be inferred from equation (4.28) of reference [43]. They read 

/N(L;0, «,<*,<*) = K + 5f L (c 1 ,c 2 )] (d - 1)/2 (3.62) 

and 

O.t^.oa) = 1 ^^ A l-i K + ^l(c 1; c 2 )] ^ 3 A 1)ljl)1 (c 1 L, c 2 L) (3.63) 
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with 

A d =^- d = -{An)~ d ' 2 r(l - d/2). (3.64) 

Building on the reorganisation of perturbation theory just described, we now wish 
to compute the residual free-energy density f TCS (L; 0,u, c±, c 2 ) and express it in terms of 
renormalized variables to obtain improved results for its scaling function D(ci,c 2 ) and 
that of the Casimir force. Clearly, from these results all Taylor expansions in e reported 
in sections 3.2 and 3.4 must be recovered, namely, the e expansions (3.48)-(3.49) of 
the scaling function D(ci,c 2 ) when min(ci + c 2 ) > 0, as well as those of the Casimir 
amplitudes A'- 0,0 -' and A^°' sp ^ given in equations (3.21) and (3.22), respectively. In 
addition, it must yield the fractional e expansion (3.25) of A^ sp ' sp ^ to 0(e 3 / 2 ) since it 
reduces to the scheme used in references [46] and [43] when c\ = c 2 = r = 0. 

What makes this approach technically difficult to implement is that the eigenvalues 
k^, whose behaviours for small Ci+c 2 we gave in equation (2.24), are not explicitly known 
for general values of b\ and c 2 . Hence we shall have to resort again to numerical means. 

Let us start by taking a look at the effective two-point vertex 



(2)/ i\z 8 2 n e f[[(p} 
r >(y-y)$af; 



(3.65) 

<p=0 



8ip a {y)5(fp{y') 

whose one-loop approximation appears in the expressions (3.62) and (3.63) for and 
fjp\ Its renormalized counterpart is given by 7^ = Z^ 2 \ Thus its Fourier p-transform 
satisfies the RG (V^ — r])^, by analogy with equation (2.39). Solving this at p = 
and r = yields the scaling form 

= 0, L, r = 0, u, Cl , c 2 , fi) » ^[El(u)] 2 L- 2+r > Q( Cl , c 2 ), (3.66) 

where Cj denote the scaling variables defined in equation (3.16). Further, is a 
non-universal amplitude whose representation as a trajectory integral can be found 
in equation (3.85c) of reference [22] but will not be needed in the remainder. Since this 
amplitude drops out of the Casimir force, we need not keep track of it and may therefore 
set it to 1 henceforth. 

(2) 

In Appendix C we determine 7^ to one- loop order, show that it is UV finite and 
compute the scaling function f2(ci,c 2 ) to 0(e). The result is 



Q(ci, c 2 ) = k 1 + eXiKi 



n + 2 
n + 8 



2 ^.2cr-l r>{a,p) TP 2 



r>\v,V to * „ i „ 

E K 1 fcig ^ Cl ,c 2 CjlllCj 7T 

(r 2 4- JVr? 4- K *\ +Kl 2^ 



, {C 2 + K 2 ){C 2 + K 2 ) L ^C 2 + K 2 2 / x 

-a,p=l V 1 1/V 2 11 3 = l 3 1 (3.67) 

, 27F a / ^ n oi o^ (ci + c 2 )(c!C 2 + /t 2 ) 2 
+ ^Am.mIc c 2 ) - (1 -21n2)« l 2(c! + K!)(c| + Kf) + 0(e ), 

where «i and Ai represent the eigenvalue «i(ci,c 2 ) and the normalization factor 
Ai(ci,c 2 ), respectively. 

The extrapolation one obtains for fi(ci,c 2 ;d = 3,n = 1) by evaluating the above 
O(e) result at e = 1 is depicted in figure 2. Since the eigenvalues «i(ci,c 2 ) are 
analytically known for all combinations of (ci,c 2 ) with Cj = 0, oo, j = 1,2, the e 
expansions of the corresponding limiting values of Q can be determined exactly. The 
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ci 1.0 
ci+1 



Figure 2. Scaling function f2(ci,C2) for n = 1 and d = 3, obtained by extrapolation 
of the 0(e) result (3.67). To cover the full range (0, oo) of the scaling variables Cj, we 
plotted Cj/ (1 + Cj) along two axes. 



same applies to the term linear in Ci + c 2 of the Taylor series expansion of f2 in Cj at 
(01,02) = (0,0). One obtains 



r^i \ i 137T 2 n + 2 2 

fi(oo, oo) = tt 2 + e— — + e 2 , 

6 n + 8 

n(o, oo) = n(oo, o) = T + + 0(e 2 ) 



(3.68) 
(3.69) 



and 



fi(ci, C 2 ) = Ci + C 2 + 0(c l5 C 2 , CiCi 



+ e 



n + 2 f 7r 



n + 8 6 



+ (Ci + c 2 ) 



7e + 1h(2tt) 



(3.70) 



+ 0(ct,c^,cic 2 ) } +0(e 



The first term on the right-hand side reflects the small-Cj behaviour k 2 (ci,c 2 ) ~ 
Ci + c 2 implied by equation (2.24). The remaining terms originate from the expansion 
of the shift 8f to linear order, 



St l (ci, c 2 ) = 5f s ^ + (ci + c 2 ) 5r£ sp + 0(c 2 , c 2 , cic 2 ), 
whose zeroth-order term is 

= ^(0,0) = ^^ r( 3-rf/2) ' 

according to equation (14) of reference [46] (where it was denoted as 5fspl P ). Note 
that this result can easily be recovered from equation (3.57) using the fact that 



(3.71) 



(3.72) 
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«m(0, 0) — (m— l)n and Ai ) i jmjm (0, 0) = 1 + 5 m ,i/2. The coefficient 5t£ sp is computed 
in Appendix D. It reads 

^, sp _ (n + 2)£iV d (d - 2)(2d - 7) T(3/2 - d/2) 

6t l ~ 3jLd _3 47r»/=»--r(3-d/2) U j - ( j 

The above results, equations (3.71)-(3.73), are consistent with equation (3.67). To 
see this note that Sr£ sp has the Laurent expansion 



L 



3L 



l + 6( 7£ -^ + ln^+0( e 2 ) 



(3.74) 



(2) 

Just as in our calculation of 7^ , the pole term gets cancelled by the counter-term 
provided by the contribution linear in b\ + c 2 = iiZ c {c\ + c 2 ) to k\. If we substitute 
equation (3.71) into k\ + 8tl along with equations (3.72) and (3.73), express the result 
in terms of renormalized variables and set u = u*, we recover the expansion of Q(ci, c 2 ) 
to linear order in Cj given on the right-hand-side of equation (3.70). 

We now insert the scaling form (3.66) of the effective two-point vertex into 
the expressions (3.62) and (3.63) for and Expressing the sum in terms of 

renormalized variables yields 

U,r( L ' °> M *> c l> C 2> A 1 ) = ^( c l> c 2) 

Th 

= -^[fi(ci,c 2 )]^ 2 (3.75) 

+ U ~^wr ^- l[n(ci ' c 2 )] d_3A M,i,i(ci, <*) + .... 

Here the ellipsis stands for contributions that are 0(w*e, e 2 ) as long as fi(ci, c 2 ) = 0(e°). 
This condition is fulfilled whenever ci + c 2 > 0. Depending on whether this is the case 
or not, the function D v (ci, c 2 ) has an expansion in integer powers of e, namely 

.3 , ^2 , 



Ap(ci,c 2 ) = - -f- + e 



12tt 



^(37 E -8 + 3.,^) 



n + 2 k 2 /ti Qq 

n + 8 8 87r 



(3.76) 

+ 0(e 2 ) ( Cl + c 2 >0), 



while 



D v (0, 0) = a 3/2 (n) e 3 / 2 + o(e 3 / 2 ). (3.77) 

Here Q' Q = 9fi(ci, c 2 ; e)/9e| e= i and a.3/ 2 (n) are the coefficient of the O(e) term in 
equation (3.67) and the quantity given in equation (3.26), respectively. 

One caveat should be noted: Although the right-hand side of equation (3.75) is UV 
finite for e < 1, it is not so at d = 3 since A d _i = [2ir(d — 3)] _1 + 0(1). The origin of 
these UV singularities is obvious. The free-energy density f v of the effective (d — 1)- 
dimensional theory at d = 3 involves UV singularities of the form A 2 and filnA. The 
bulk and surface counter-terms included so far do not cure these; additional subtractions 
(of the kind needed for a super-renormalizable effective bulk theory in d — 1 = 2 



24 



dimensions, cf. [48]) would be needed. Evaluated at d = 3, the first term on the right- 
hand side of equation (3.75) would then yield a contribution (Q/8-7r)[l — 7s — ln(f2/47r)] 
to / v ,res,R and hence to D v . Such logarithmic terms are encountered at d = 3 also in 
Gaussian film models; see e.g. reference [49]. 

Let us postpone any further discussion of the behaviour of D v at d = 3 for the 
moment and first compute the free-energy contributions and fT 1 . The former differs 
from /M through the lowest-mode contribution. This is given by equation (3.62) with 
the shift 5tl set to zero. Hence we have 



4 11 (L;0,w > ci,c !2 )-/L 1 J(L;0,«,c 1 ,c 2 ) 



,d-l 



d 



(3.78) 



[21 

To determine the two- loop graph of , we must subtract from the two- loop graph of 
/™ displayed in equation (3.29) the contributions for which the mode indices m x and 
m 2 of the upper or lower loops take the value 1, i.e. 
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Details of the calculation are presented in Appendix E. They yield 

L d-i 



(3.79) 
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(3.81) 



A/,,i( c i! c a) = A(ci, c 2 ) + —4 (8 - 3 7E - 3 ln(«iA0) 
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(3.82) 



In 4 (ci + c 2 )(cic 2 + k?) 



2 - ( C ? + K?)(cf + K?) 

To visualise the result, we have plotted in figure 3 the difference between the 0(e) 
expression (3.48) for D(ci,c 2 ) and its analogue (3.80) for D^,(ci,c 2 ) one obtains in 
the scalar case n = 1 upon evaluation at e = 1. This difference corresponds to the 
contribution from the lowest mode m — 1 to the series expansion of D(ci,c 2 ) to 0(e). 
It vanishes as (ci, c 2 ) approaches the origin and reaches its minimum at the fixed point 
(ci,c 2 ) = (00,00) corresponding to large-scale Dirichlet boundary conditions at both 
boundary layers 53 1 and 53 2 . 

Combining equations (3.75) and (3.80), we can write the result of the modified 
RG-improved perturbation theory used in this subsection as 



D(ci, c 2 ) = D v (ci, c 2 ) + D^(ci, c 2 ) 



(3.83) 
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Figure 3. Contribution from the lowest mode m = 1 to the e expansion of the 
scaling function _D(ci,C2) for n = 1, evaluated at e = 1. This quantity is given by 
{D (c 1 ,c 2 ) - A/,,o(ci,c 2 ) + e[L>i(ci,c 2 ) - A/,,i(ci, c 2 )]} e =i,n=i. 

where stands for the expression defined through the right-hand side of equation (3.75) 
in conjunction with equation (3.67), while represents the 0(e) series expansion in 
equation (3.80) with the expansion coefficients given by equations (3.81) and (3.82), 
respectively. Note that the result must be utilised with care. It is not suitable for direct 
evaluation at d = 3 (e = 1). We defer a discussion of this and related issues to section 4. 

3.5.2. Consistency with fractional e expansion for C\ = c 2 = and Ginzburg-Levanyuk 
criterion Our motivation for working out the modified RG-improved perturbation 
theory described in section 3.5 was the goal to achieve consistency of the theory for 
general nonnegative values of the enhancement variables c\ and c 2 with the approach 
used in references [46] and [43] to study the zero-mode case c\ = c 2 = of critical surface 
enhancements. Since the former approach reduces to the latter when c\ = c 2 , it is clear 
that consistency is ensured and the fractional e expansion (3.25) must be recovered. To 
see this explicitly from the results of the previous subsection, recall that Ki vanishes 
as Ci + c 2 — > 0. Therefore, D^ !0 (ci,c 2 ) and D^ i i(ci,c 2 ) approach the limiting values 
D (0,0) and -Di(0, 0), respectively. Furthermore, the first term on the right-hand side 
of equation (3.75) yields the e 3//2 term according to equation (3.77). 

That the modified RG-improved perturbation theory is controlled for small e, 
irrespective of whether c\ + c 2 is positive or zero, can also be seen from a criterion of the 
Ginzburg type [50, 51]. Let 7^(pi,P2,P3) {^ 1T ) d l ^Cl2t=i Pi) ^ e ^ ne Fourier transform 
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of the effective four-point vertex 7 ( - 4 * ) (j/i, . . . , 2/4). The analogue of equation (3.66) tells 
us that (0,0,0) ~ u*fi 2ri L d 5+2, ? ) where the proportionality factor is a function of 
the scaling variables Ci and c 2 . Following a standard reasoning (used also by Sachdev 
in his study of finite-temperature crossovers near quantum critical points [48]), we 
can construct from 7^(0), its derivative d p 2^ 2 ^\ and 7^(0,0,0) the dimensionless 
interaction constant 

f 4 ) (0,0,0) 

U = [(<9 p27 (2) ) P =o] (d - 1)/2 [7(2) (0)] (5-^0/2 • (3 ' 84) 

Since T 2 '(0) is proportional to Jl(ci,c 2 ), this measure of the strength of the 
nonlinearities in "H e ff behaves for small e as 

r, * ( e (1 ~ e)/2 if ci = c 2 = 0, 
[ e if ci + c 2 > 0. 

Hence the Ginzburg-Levanyuk-type criterion U <C 1 is satisfied when e < 1. 



4. Discussion of results and extrapolation to d = 3 dimensions 

As our analysis in section 3 has explicitly confirmed, the scaling functions D(ci,c 2 ) 
and T>(ci, c 2 ) of the residual free-energy density and the Casimir force have power series 
expansions in e when Ci+c 2 > 0. Their coefficients of the corresponding series expansions 
to first order in e are given by equations (3.45)-(3.49) and follow from equation (3.18), 
respectively. Assuming that Ci + c 2 > 0, the most elementary extrapolation to d — 3 
dimensions we can make is to set e = 1 in these O(e) expressions. In figure 4 a contour 
plot of the resulting extrapolation for the scaling function D(c 1 ,c 2 ) is depicted as a 
function of the variables Cj/(cj + 1) with j = 1,2. A three-dimensional plot of this 
result can be found in figure 1 of reference [21]. The plotted function D is negative 
in the vicinity of the 11 diagonal and changes sign across the thick dashed lines. The 
blue lines indicate paths generated by changes of L at fixed values of the enhancement 
variables C\ and c 2 , i.e. RG flow lines. Note that the path associated with a given pair 
(ci,c 2 ) of enhancement variables does not depend on the value of the surface critical 
exponent y c > appearing in the definition of the scaling variables Cj in equation (3.16) 
once the non-universal amplitude E*(u) has been fixed. As is borne out by figure 4, 
choices of (ci, c 2 ) exist for which these paths start in the region D < 0, enter the region 
of positive values as L increases and finally return back to the region D < 0. Hence 
one expects that crossovers of the Casimir force from attractive to repulsive and back 
to attractive behaviours should occur as a function of L. 

The result for the scaling function X> n= i(ci, c 2 ) of the Casimir force that follows from 
this extrapolation for D n=1 (ci, c 2 ) via equation (3.18) is displayed as a three-dimensional 
plot in figure 5 and as contour plot in figure 6. 

To obtain these plots, we proceeded as follows. We substituted the extrapolated 
expression [D + eDi] €= i for the scaling function D in equation (3.18). For the pre- 
factor d — 1 + y c we used the value 2.718 corresponding to Hasenbusch's recent Monte 
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C2 
C2 + 1 




Cl+1 

Figure 4. Contour plot of the 0(e) result (3.48) for the scaling function D n= i(ci, C2) 
of the residual free-energy density / ros at the bulk critical point, evaluated at e = 1. 
Along the axes the variables Xj — Cj/ (Cj + 1) are plotted. For further explanations see 
main text. 




Figure 5. Three-dimensional plot of the 0(e) result for the scaling function T> n= i 
[equation (3.18)] of the Casimir force at T C OD , evaluated at e = 1. The axes are chosen 
as in figure 4. The thick dashed lines mark the locations of zeros of T>. 
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ci+1 



Figure 6. Contour plot of the 0(e) result for the scaling function T> n= i 
[equation (3.18)] of the Casimir force at T Ci00 , evaluated at e = 1. The axes are 
chosen as in figures 4. The thick lines mark the locations of zeros of T>. The blue lines 
indicate the flow induced by changes of L at fixed surface enhancements (ci, ci) 

Carlo estimate y c (n=l, d=3) ~ 0.718(2) [52]. Let us note that some earlier Monte Carlo 
calculations led to significantly larger values for this exponent [53, 54, 55], but others and 
more recent ones [56, 57, 58, 59] yielded similar numbers. Field-theory estimates based 
on the e expansion to second order and the massive field-theory approach to two-loop 
order [60, 61] gave y c (l,3) ~ 0.75 and y c (l,3) ~ 0.85, respectively. For a detailed list of 
Monte Carlo and field-theory estimates for y c (l, 3) the reader is referred to reference [52]. 
Choosing a somewhat larger value of y c (l, 3) such as the quoted field-theory estimates 
would lead to small quantitative, but no qualitative changes in figures 5 and 6. 

The behaviour of the plotted function T> is qualitatively similar to that of D. As 
is obvious from figures 5 and 6, the critical Casimir force does exhibit the anticipated 
crossovers from attraction to repulsion and back to attraction for appropriate choices 
of (01,02). According to a theorem for systems satisfying reflection positivity [62], the 
Casimir force cannot become repulsive for equal enhancements c\ = C2- This conforms 
with the fact that the Casimir force, by continuity, is attractive (i.e. T> < 0) in the 
vicinity of the 11 diagonal. By choosing (01,02) sufficiently away from the 11 diagonal, 
one can ensure that a double sign change of J 7 occurs as L increases. A fascinating 
consequence is that for such choices of (01,02), critical values L = L (ci,C2) of the film 
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thickness exist for which the critical Casimir force vanishes. 

In the extrapolations for D n=1 (ci,c 2 ) and £> n=1 (ci,c 2 ) presented above, we tacitly 
assumed that c\ + c 2 > 0. Knowing that the e expansion breaks down when c\ = c 2 = 
(becoming fractional), we may expect to see indications of this breakdown in the 
behaviour of these extrapolations near the origin. This is indeed the case, though only 
very close to it. To demonstrate this, we depict in figure 7 the values [D + eDi] n= i )e= i 
of the naively extrapolated scaling function (3.48) on the diagonal Ci = c 2 = c. 




Figure 7. Extrapolated 0(e) results [Do + A]n=i(c, c) (red curve) and [D^fl + 
D^ t i] n =i{c, c) (blue curve); cf. equations (3.48) and (3.80). The brown line is the 
function D app specified in equation (4.2). The asymptotic small-c form (4.1) of the red 
curve is plotted as dashed red line in the in-set. The dotted horizontal lines indicate 
the limiting coinciding values of the two red and blue curves, and that of the brown 
curve, respectively. For further explanation, see main text. 



This function [D + Z?i] n= i(c, c) has an infinite slope at c = 0. Its asymptotic 
behavior near the origin can be determined analytically in a straightforward fashion. 
One finds 

r 2 



[D + D 



7T 



1 71=1 



c, c 



+ ai(l 



7T 

144 



V^c + Ofc) 



>o 1440 

: -0.0116593 - 0.0308534 y/c + O(c). 



(4.1) 



The in-set in figure 7 shows a comparison of [Do + -Di] n =i(c, c) and its asymptotic 
form (4.1) for small c. The term oc c 1//2 in equation (4.1) results from that part of the 
two-loop contribution flel whose mode summation ^2 mi m2 is restricted to vn\ = 1 in 
either one of the two loops and to m 2 > 1 in the respective other. The easiest way to 



recover it is to expand D„ 



c,c 



-_i to O(e). The coefficient of the O(e) term is precisely 



the contribution oc a/c in equation (4.1). It is a consequence of the (spurious) infrared 
singularity one encounters in this expansion due to the m = 1 mode with K\(c, c) — > 0, 
and would imply an infrared singular derivative d c of the function (4.1) at c = 0. This 
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is a spurious infrared singularity since the system is not expected to exhibit critical 
behaviour at T — T c oo = c x = c 2 = when L < oo. The extrapolated scaling function 
[Do + ePi] n =i je= i of the Casimir force behaves analogously near the origin, as should be 
clear from equation (3.18). 

To illustrate that the spurious behaviour of [D + Di\ n=1 (c, c) near the origin is 
due to the zero-mode contribution, we also display the function [A/,,o + A/>,i]n=i(c, c) 
in figure 7. The latter function is regular at c = and has a finite positive slope there, 
since it does not involve the zero mode m\ — \. Its value at c = agrees with that of 

[A) + £>l]n=l- 

The remaining brown curve in figure 7 represents an extrapolation based on 
equation (3.83), namely the function 

D app (c, c) ~ (lV, (c, c) + eD^c, c) - (A 3 /3) [Q(c, c; e)] 3/2 | , (4.2) 

l. J n=l,e=l 

This means that we have discarded the two-loop contribution to D v and replaced the 
pre-factor Ad-i/(d — 1) of the one-loop term by its value at d — 4. To understand the 
rationale for this approximation, one should note the following. If Ci = c 2 = 0, then 
this term reduces precisely to the a3/ 2 (n)e 3 / 2 contribution to D(0,0). If we expand, on 
the other hand, this choice of D v to linear order in e when ci + c 2 > and substitute it 
into D^ o + eD^ x + D v , we recover the 0(e) result D + eD 1 . 

Hence, in the limit c — > 0, D app (c, c) approaches the value one obtains for 
A( sp ' sp ) | n= i j£= i by evaluating the small-e expansion (3.25) at e = 1. Furthermore, its 
derivative d c D app is finite at c = 0. Despite these appealing features, we refrain from 
showing a plot of this function for all c G (0, oo) because away from the origin this 
function deviates considerably from the O(e) extrapolation [D + D]\ n=1 . The reason is 
that f2(ci, c 2 ; e = 1) becomes fairly large for large Ci +c 2 (see figure 2) and the difference 
between (A 3 /3)f2(c, c; e) and its extrapolated O(e) expansion (l+ed e )(A 3 /3)Q(c, c; 0) | e=i 
is not small. Away from the origin, we therefore do not consider D app to be superior to 
the naive extrapolation D + D\. 

As we mentioned earlier, the pre-factor A&-\j id — 1), has a UV pole at d — 3. Being 
guided by the aim to achieve consistency with the small-e expansion for ci = c 2 = 0, 
we expanded it in e and hence could replace it by its value at d — 4 at the order of 
our calculation. In an extrapolation based on equation (3.83) one might be tempted 
to use the value of the UV finite one-loop term of D v one obtains in the manner 
described below equation (3.77) by subtracting the pole term oc (d — 3)" 1 from 
— ^3j-[Q(ci, c 2 , e = l)]*^ 1 )/ 2 and evaluating the difference at d — 3. However, such 
an approach would mix the RG approach based on the e expansion we used for with 
elements of a fixed-dimension RG approach. The corresponding extrapolation would 
again lead to substantial differences from the O(e) extrapolation D + D\ for large 
values of C\ + c 2 . We feel that such attempts to evaluate D v directly at d = 3 should 
be based on a RG approach in fixed dimension d = 3. Such an approach is known to 
involve the study of the theory away from the bulk critical point [63, 64]. Furthermore, 
the surface enhancement variables c\ and c 2 provide two mass parameters in addition to 
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the bulk correlation length one has to deal with. A corresponding massive RG approach 
in fixed dimension d = 3 based on references [60] and [61] — or combined with elements 
of the strategy followed in reference [65] - is technically very demanding and beyond 
the scope of the present work, albeit conceivable. 

It is instructive to compare the extrapolated scaling function V n=1 of the critical 
Casimir force displayed in figures 5 and 6 with their exact analogues 

ZWci, c 2 ; d) = (d - 1 + ci9 Cl + c 2 9 C2 )/W(l; 0, ci, c 2 )/n (4.3) 

of the Gaussian theory in d — 3 and d = 4 dimension. Here, the one-loop term (3.41) 
must be substituted for fill. A plot of these functions on the diagonal Ci = c 2 = c is 
shown in figure 8. 
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Figure 8. Comparison of the scaling functions £>Gaufi(c, c; d) with d — 3 (dashed) and 
d = 4 (dash-dotted) from equation (4.3) with the extrapolated function V n=1 {c,c) 
displayed in figures 5 and 6 (full line). 

As one sees, the extrapolation T> n= \ (c, c) lies for most values of c below and above 
the Gaussian scaling functions (4.3) with d = 4 and d — 3, respectively. Having 
already discussed the spurious small-c behaviour of the extrapolation, let us add a 
comment about the limiting behaviour at large c. It has been known for long that 
deviations of the surface enhancement variables Cj from their fixed point values Cj = oo 
associated with the ordinary surface transition correspond to irrelevant surface scaling 
fields ("extrapolation lengths") whose RG eigenvalues are given by the momentum 
dimension —1 [66, 22]. Taking into account the corrections to scaling oc L~ x implied 
by them is crucial for a proper analysis of Monte Carlo simulation results for Casimir 
forces [26, 27, 28, 29, 30]. Our small-e results enable us to explicitly verify the presence 
of such corrections to scaling. In order to comply with a correction to scaling to the 
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residual free energy that is down by a factor 1/L, the scaling function D(ci,C2) ought 
to exhibit the asymptotic behaviours 

, (0,0) 



C1,C2— >oo \ 
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Inspection of our results given in equations (3.45)-(3.49) reveals that contributions 
~ eln(Cj)/Cj originate from Di which can be exponentiated in accordance with these 
limiting forms. The universal amplitudes A^' ^ and /\f p '°^ are found to have the e 
expansions 
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respectively. Let us recall that the behaviour of the extrapolation T> n= \ plotted in 
figure 8 was obtained by substituting [D + eDi] e=1 into equation (3.18). Therefore, 
the plotted extrapolated function does not approach its limiting value at c = 00 as 
c -1 /^. This is because its asymptotic form contains the logarithmic term ~ em(c)/c 
mentioned above, evaluated at e = 1. Hence, the deviation of the plotted extrapolation 
V n= i from its value at c = 00 varies (incorrectly) as ~ ln(c)/c. For conciseness, we have 
refrained from designing alternative extrapolations into which the correct asymptotic 
large-c asymptotics is incorporated via appropriate exponentiation ansatzes. 



5. Summary and Conclusions 

Considering an 0(n)-symmetric 4 model in film geometry, we investigated the 
fluctuation-induced forces produced by thermal fluctuations between the boundary 
planes of a film that is held at its bulk critical point. We restricted our attention to the 
case of generic non-symmetry-breaking boundary conditions. The interaction constants 
&j of the 0{n) symmetric boundary terms Cj0 2 /2 included in the Hamiltonian were 
chosen to correspond to non-supercritical surface enhancements but otherwise allowed 
to take arbitrary values. Together with the restriction to the bulk critical temperature 
T c> oo (or, more generally, to temperatures T > T c oo ), this requirement guarantees that 
the 0(n) symmetry is neither explicitly nor spontaneously broken. 

Whenever the symmetry <fi — > —<fi is broken, the order-parameter profile (0) does 
not vanish. As discussed in the introduction, this implies that a nonzero Casimir force is 
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obtained even in mean- field theory — i.e., in the absence of fluctuations. Thus in these 
cases the Casimir force is not fully fluctuation induced but contains a non-fluctuating 
component. This applies, in particular, to binary liquid mixtures, for which beautiful 
direct measurements of the Casimir forces were accomplished recently [16, 29, 67]. A key 
motivation for excluding such symmetry breakdowns in the present investigation was 
the intention to keep the thermodynamic Casimir forces entirely fluctuation induced, 
making them share this crucial feature with the original QED Casimir force [1, 4]. 

Our analysis was based on the field-theoretic RG approach in 4 — e dimensions. 
The use of the e expansion in studies of critical and near-critical Casimir forces has 
shortcomings that go beyond the familiar ones known from its application to critical 
behaviour in bulk and semi-infinite systems. The series expansions it yields for bulk and 
surface critical exponents as well as for other properties at bulk critical points — Taylor 
expansions in e — are asymptotic. Naive extrapolations of such expansions to low orders 
in e are not normally quantitatively reliable. However, quantitatively accurate results 
can be obtained by extending the expansions to sufficiently high orders, combining them 
with information from appropriate large-order results and using appropriate Borel-Pade 
resummation techniques [68, 64, 69, 70]. The additional complication one encounters 
in the case of films is that e = 4 — d ceases to be an appropriate expansion parameter 
when dealing with dimensional crossovers. (The adequate analogue of this parameter 
for the study of critical behaviour at the transition temperature T Cj l of a film of finite 
thickness L and infinite extension in d— 1 Cartesian directions would be e 5 = 5 — d.) The 
limitations of the approach become particularly apparent when zero modes appear at 
zero-loop order at the bulk critical point, as happens for critical enhancement c\ = C2 = 
of both boundary planes. Whenever this occurs, the e expansion breaks down at the 
bulk critical temperature T Cj00 and becomes fractional, involving also powers of lne 
[43, 46, 47, 48]. By contrast, the e expansion remains (formally) intact at T c>00 provided 
at least one of the (non-supercritical) enhancements Cj is subcritical, although a similar 
breakdown clearly must occur at a temperature below T Cj00 . 

Despite these deficiencies, the field-theoretic RG approach in 4 — e dimensions has 
a number of appealing and valuable features. First of all, it lends itself to nontrivial 
checks of the scaling forms the residual free energy and the Casimir force are expected 
to have according to finite-size scaling theory and the RG approach. Since we performed 
a two-loop calculation which gave / rcs to order u, we were able to identify the ln(/xL) 
contributions implied by the nontrivial 0(e) term of the surface crossover exponent 
<3> in the scaling arguments (3.13) and thus verify the scaling forms (3.13) and (3.17) 
to first order in e. Second, the fact that general properties such as the exponential 
decay in the disordered phase (see references [20, Sec. VI] and [43, Sec. IV. C]) and 
analyticity properties of the free-energy density / (see references [20, Sec. VII] and 
[43, Sec. IV. C]) must hold at any order of the small-e expansion provides nontrivial 
checks of them. For the sake of clarity, it will be helpful to recall these analyticity 
properties. For given n and sufficiently large d, the film undergoes a sharp transition to 
an ordered phase at a temperature T c ^ < T C)OC . In the case of non-supercritical surface 
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enhancements considered here, there should also be no surface transitions of the film 
at temperatures T > T coo (i.e. its surface layers should not exhibit long-range order). 
Consequently, the total free-energy density / must be analytic in temperature at T Cj00 . 
This requirement imposes conditions on the behaviour of the temperature dependent 
analogue of the scaling function of / res . We refer the reader to references [20, Sec. VII] 
and [43, Sec. IV.C] for details. 

These analyticity requirements were found to be fulfilled by the two-loop e- 
expansion results for anti-periodic, (O, O), and (O, sp) boundary conditions, but violated 
at 0(e) for periodic and (sp, sp) boundary conditions [20]. Extension of the small-e 
expansion to 0(e 3 / 2 ) in the latter two cases via the effective-action approach outlined 
above increased the order at which a violation of analyticity requirements occurred to 
e 3/2 |43j_ Third, our analyses in the latter reference and above yield evidence of a shift 
of the point c\ = = at which the zero-mode field tp becomes critical at T Cj00 to 
an L-dependent location with Ci + c 2 < 0. By analogy with the shift T c oo — > T c L of 
the transition temperature, such a shift is expected and in conformity with the required 
analyticity of / at bulk criticality. 

One crucial problem one is faced with in studies of critical behaviour and 
dimensional crossover in films is the difficulty of obtaining accurate results for such shifts 
in d = 3 bulk dimensions by analytical methods. The reason is that fluctuations on all 
length scales between the microscopic one (lattice constant) and the bulk correlation 
length may contribute, as a result of which these shifts may acquire contributions that 
are non-analytic in the interaction constant [71, 72, 60, 73]. Since we used RG-improved 
perturbation theory in 4 — e dimension at T Cj00 to determine the effective action T-L e s, our 
approach is not capable to capture such non-analytic contributions. It gave us a shift 
along the ci+c 2 axis of the point where the ip component of the order parameter becomes 
critical that was proportional to u*. This in turn implied a contribution to the critical 
Casimir force for critical surface enhancements c\ = c 2 = of the form (u*)^ 3 ~ £ ^ 2 . To 
make our results for the scaling functions D(ci,C2) and X>(ci,C2) consistent with the 
fractional e expansions at (ci,C2) = (0,0), we extended the effective-action approach 
of references [46] and [43] to nonzero values of ci + C2. We succeeded in achieving 
consistency with the results of these references for the case c x = c 2 = and the fractional 
e expansions. Furthermore, the effective-action approach turned out to be well-behaved 
for small e. Nevertheless, its results must be taken with a grain of salt: they did not lend 
themselves easily to simple extrapolations to d — 3 dimensions that could be judged as 
clearly superior to those based on the conventional e expansion. 

Let us emphasise that our results exhibit a number of interesting and important 
qualitative features which may be expected to persist in quantitatively more accurate 
investigations, irrespective of the accuracy of our extrapolation to d — 3. 

(i) Generalising previous work [19, 20, 37, 43, 47], they clearly demonstrate that critical 
Casimir forces can be attractive or repulsive depending on the values of the surface 
enhancement variables Cj, even when the internal symmetry (Z 2 for n = 1 and 0(n) 
for n > 1) is neither spontaneously nor explicitly broken. 
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(ii) They are in full accordance with and hence corroborate - the scaling 
forms (3.14) and (3.17) of the residual free energy and Casimir force, respectively. 
This means that the Casimir amplitude becomes scale dependent. 

(iii) As a further important observation we found that crossovers from attractive to 
repulsive forces and vice versa can occur if the surface interaction constants take 
appropriate values. In conjunction with (ii) this means that the critical Casimir 
force goes through zero at a certain thickness Lq. 

(iv) For the case of primary interest — the Ising bulk universality class with short-range 
interactions and d = 3 — we saw that for properly chosen surface enhancement 
values Ci and c 2 , crossovers from attractive to repulsive and back to attractive 
behaviours should occur as L increases. It would certainly be worthwhile to check 
this prediction along with (iii) by Monte Carlo calculations. 

(v) Finally, let us mention that crossover behaviours analogous to (iii) and (iv) should 
also be possible for three-dimensional 0(n) systems with easy- axis spin anisotropies. 
Semi-infinite systems of this kind are known to have anisotropic analogues of the 
isotropic special transition [74, 75], characterised by the coincidence of the transition 
temperature at which the easy-axis component of the order parameter at the surface 
becomes critical with the transition temperature T c oo of the n > 1 bulk system. 
Hence for finite thicknesses L, the surface interaction constants associated with 
a given easy axis on one or both boundary planes can be critically enhanced. 
This suggests that the mentioned analogues of (iii) and (iv) should occur. (Films 
involving different directions of easy axes on the two boundary planes would require 
separate analyses.) 

Appendix A. Calculation of the one-loop free-energy density /M 

Appendix A.l. Calculation based on the Abel-Plana summation formula 

To compute /W, we start from equation (3.28), differentiate this expression with respect 
to r and compute the integral over p. This gives 



9 f /W(L;f,c 1 ,c 2 ) = ^ / 

1 m=l J P 



p T + P 2 + 

(A.l) 



m 

oo 



= _ L -{d-3) rmKd-i Y"(tL 2 + k 2 ) {d ~ 3 ^ 2 
4cos(d7r/2) l^ TL +Km) 

The expression in the second line of equation (A.l) varies ~ L J °° dk (f + &; 2 )( d-3 )/ 2 
in the large- L limit. We could subtract this bulk term to make the difference well 
defined. However, we are ultimately interested in the free-energy density f^(L; f, c±, c 2 ). 
Term-by-term integration of the series for increases the exponent (d — 3)/2 of the 

series coefficients by one. Hence, subtracting the limiting bulk contribution Lf^ would 
not suffice to render UV finite expressions for the free-energy density f^(L; f, ci, c 2 ) 
in dimensions d < 4. To obtain well-defined results for this quantity, we prefer to use 
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analytical continuation in d, rather than making additional subtractions. We start by 
integrating the series in equation (A.l) term-wise, choosing the integration constant 
such that 

^ *■ - « - g (fL2 + K ™ )< "" /2 <a ^ 2) 

That the integration constant has been chosen correctly may not be obvious at 
this stage. However, once we will have analytically continued the above series (A. 2) 
in d, one can confirm this choice by verifying that the resulting free-energy expressions 
reduce to known dimensionally regularised results at f = 0+ in the limits 61,62 — > 
and 61,62 —> 00. Furthermore, there are two other ways to arrive at expression (A. 2), 
both of which indicate its consistency with general rules of dimensional regularisation. 
The first is to use the representation 

\ nv = ™ (e-*-e-^) (A.3) 
Jo t 

for the logarithm in the second line of equation (3.28), perform the integrations and 
take into account that integrals of pure powers such as ^ 1 vanish in dimensional 
regularisation. Another way is to insert the partial-p identity (d — 1) _1 V P • p into the 
integrand of the integral f^ d ^ in the last line of equation (3.28) and integrate by parts, 
dropping the boundary terms at p — 00. 

According to the result (A. 2) we must calculate — and analytically continue — a 
series of the form 

00 

S ClA M) = 5>^ + 6 2 r (A.4) 

m=l 

with b 2 = fL 2 , where the positive solutions to equation (2.21). Such series are 

generalisations of Epstein-Hurwitz ( functions (see e.g. reference [76]). Let us assume 
that both Cj > and C\ + C 2 7^ 0. Then all K m , m = 1, 2, . . . , 00, are positive, and as 
we showed in section 2, non-degenerate. Hence the function (9 K lni?Ci,c 2 ( K ) nas a simple 
pole with residue 1 at each of these n m . Since (k 2 + b 2 ) a is regular at k — n m , the value 
of this function at n m is given by the residue Res \(k 2 + b 2 ) a B! c C2 {k) / Rc^i 1 ^)]- But 
at k = K m R' c Ca can be factorized as 

R'c^cMm) = 2A r Cl ,c 2 ( K m)/^ = N Cu c 2 (K m ) g Cl ,c 2 (K m ), (A.5) 



where 



while 



JV CiA (k) = (dC 2 - k j cosk- (C 1 + C 2 )K,smK, (A.6) 



9c u c 2 (k) = 1 + ^— — + -2— — . (A.7) 



Ci | c 2 

Thus we have 

{K 2 m + b 2 ) a = R e s[(K 2 + b 2 rr Cl ,c 2 (K)] (A.8) 
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with 



T Cx,C 2 ( K ) 



-Rci,c 2 ( k ) 



(A.9) 



Since this function and the one appearing in Sc lt c 2 ( a ! b) are even in k, we can extend 
the summation in equation (A. 4) to all integer- valued m^O, using K-\ m \ = — K[ m |, and 
divide by 2. Since (k 2 + b 2 ) a is regular at /c m , each term of this series can be expressed 
as (27ri) _1 times an integral § dn (k, 2 + b 2 ) a Tc lt c 2 { K ) along a contour that passes once 
around the pole at the respective counter-clockwise fashion and contains no 

other singularities. If a < —1/2, the series Sc 1 fi 2 { a '^) converges and the integrand of 
the contour integral decays sufficiently fast as k — > ±oo ± iO, so that the union of all 
these contour integrals can be deformed into a path 71 encircling all poles K m with 
7^ m G Z (see figure Al). We now add and subtract integrals along the paths 72 
and 73 depicted in this figure. Since the integrand is regular at all nonzero k inside the 
region bounded by the closed path 71 U 7 2 U 73, the integral along it is — 27ri times the 
residue of the integrand at k = 0. Furthermore, the integral along 72 approaches zero 
as the radius of the circle on which 72 is located tends to infinity. Hence we have 



The factor (k 2 + b 2 ) a implies that the integrand has branch cuts from ib to ioo 
and from — ib to —ioo. In addition, the function gc u c 2 ( K ) an d hence Tc 1) c 2 («) have 
poles at ±iCi and ±iC 2 , marked by yellow pentagons in figure Al. We move the path 
73 infinitesimally close to the imaginary axes, passing around the poles ±Cj along 
semicircles of radius 6 — > 0. The portions of 73 with Imn > and Im/c < give 
identical contributions. Taking into account that the limiting values of the power on 
the right and left rim of the branch cut are 




(A.10) 




Figure Al. Paths in the complex K-plane. 




2 



.2 



2 



.2 



±i 




Re K->0± 
Im K>0 
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one arrives at 



1 f°° dt 
S Cl ,c 2 (a;b)=--b 2a -2sm(ira)V / - (t 2 - b 2 ) a T ClA (it) 

2 Jb TT 

2 

- 4th ^ (Cj - 6 2 ) a cos(a7r) Res T ClA (/c). 



(A.12) 



dt, where C< 



Here P dt means the principal value lim,5^o+ ( Jj, C< 5 + j'c<+8 Jc 
and C> are the smaller and larger one of C\ and C2, respectively, and b is assumed to 
be smaller than C<. The contribution in the second line results from the integrals along 
the semicircles. 

In order that the integral V dt exists, we must require a > —1 (to ensure 
convergence at the lower limit of integration) in addition to the original condition 
a < —1/2. In view of the behaviour of the integrand near the upper integration limit, 
we must have a < —1/2 (to guarantee convergence at the upper limit of integration). 
However, we can split off the limiting value lim^iioo Nc u c 2 ( K ) I Rci,c 2 { K ) = Th obtaining 



N Cl ,cM 

R d,C 2 ( K ) 



=Fi forlm/t^O. (A. 13) 



C\ ± in C2 ± i« 

Since the first term vanishes at k — iCj, the contribution it yields to the term in the 
second line of equation (A.12) vanishes, i.e. Res [T CljC2 (/t) + igc u c 2 ( K )] = 0- The two 

K=\Cj ' ' 

contributions to Sc lt c 2 ( a ', b) implied by the term — i in equation (A. 13) can be represented 
as an integral along the original path 7 1? namely 
1 



'CiA 



1 f f°° Hk 

a;b) = — dKsgn{\mK)g Clfi2 {K){K 2 +h 2 T= — g Clfi2 {n) (k 2 + b 2 ) a ,(A.14) 
47T J-n Jo 7r 

where sgn(x) means the sign function. For a < —1/2, the latter integral is well defined 

and can be analytically computed. This gives 



I CuC2 (a;b)=b 1+ 



2a 



1/2) 



+ 



2 v /7rr(-a) 2cos(7ra) 



e ( c J-"y 



,l+2a 



Y(-a)Cj 



2 F 1 (l/2,l;a + 3/2;b 2 /C 2 ) 



(A.15) 



where 2F1, the regularised hypergeo metric function, is an entire function related to the 
standard hypergeometric function 2 Fi via 

2 F 1 (a, 7; z) = 2 F 1 (a, /3; 7; ^)/r( 7 ). (A.16) 

The expression on the right-hand side of equation (A.15) provides the analytic 
continuation of interest of the integral (A. 14) to positive values of the exponent a. 
Our result for the series Sc lt c 2 (o; b) thus becomes 



S Cl ,c 2 (a; b)= - - b 2a + I Cl ,c 2 (a; b) - 2sin(7ra) 



dt gc u c 2 (it) (t 

71 



L2\a 



t+Ci t+C 2 2t 
t-d t-C 2 c 



(A.17) 



where Ic l fi 2 { a '- I b) represents the meromorphic function defined by the right-hand side of 
equation (A.15). 
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Equation (A. 17) provides the appropriate generalisations of the results of 
reference [20] for the special cases (Ci,C 2 ) = (0,0), (oo, oo) and (0, oo) to general 
nonnegative values of C\ and C 2 . To check the consistency with Krech and Dietrich's 
results, let us work out the asymptotic behaviours in these limits. Straightforward 
analysis yields 

dt (t 2 - b 2 ) a 



2 sin(7ra) 



7T 



3 2t 



lim 

C2-»0O 



(a 



<2\a 



2 cos(7ra) 



1/2) 



+ 2 sin(7ra) 



dt (t 2 



<2\a 



(A.18) 



(A.19) 



7T 



^2f 



+ 1 



and 



lim 

Ci,C2-*oo 



Sc u c 2 {a; b) 



U2\a 



2 cos(7ra) 



S 0>0 (a; b) - b 



2a 



(A.20) 



The right-hand sides of equations (A.20) and (A.19) correspond to equations (C5) 
and (C7) of reference [20], respectively, and our result (A.18) complies with the remark 
about the Neumann-Neumann case below its equation (C5). To see this, one should 
identify b 2 with tL 2 in their notation, multiply our result with the factor L~ 2a (= L 1_d ) 
and take into account the relation sin(7ra) = — 7r/[r(— a) T(l + a)]. The terms we have 
subtracted on the left-hand sides of equations (A.19) and (A.20) approach infinity in the 
considered limits Cj — > oo. These do not appear if one sets C 2 = oo and C\ = C 2 = oo 
from the outset because the operations of analytic continuation and of taking these 
limits do not commute. (As long as a < 0, the subtracted terms approach zero as 
Cj — > oo.) Note also that these terms contribute only to the excess surface densities /j 1 !, 
but not to fill. 

We now insert the above results (A. 15) and (A. 17) into equation (A. 2), set 
a — (d — l)/2, and take the limit b — > (i.e. f — > 0). The ^-integral of equation (A. 17) 
becomes 



g Cl ,c 2 (tt)t 



d-l 







t+Ci t+C2 2t 
t-d t-C 2 



1 



dt .d-l d 



t 







7T 



dt 



In 



(t-d)(t-C 2 )e- 
(d+t)(C 2 + t) 



and can be integrated by parts. An elementary change of variable t — > p 
leads us to the result 

n ft Kd-x 



/W(L;0,c 1 ,c 2 ) = - 



+ - 2 K d . 1 



2 (d-l) sin(7ro?) 

;>oo 

/ dpp d ~ 2 \n 
Jo 



d-i 
l 



_ (p-C!)(p-c 2 )e- 2L P 
(ci +p)(c 2 +p) 

which is in conformity with equations (3.33), (3.34) and (3.41). 



(A.21) 
t/L finally 

(A.22) 
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Appendix A. 2. Calculation based on the free propagator 

In this appendix we describe an alternative calculation of the dimensionally regularised 
one-loop free-energy density f^(L; 0, ci, c 2 ), which uses the result (2.28) for the free 
propagator. We start from the one-loop expression 

d f fW(L-T,c 1: c 2 ) = ^j { " ^ £dzG L (p-z,z\T,c 1: c 2 ) (A.23) 

for the energy density. Upon inserting the free propagator with z — z' from 
equation (2.28), the integration over z can be performed using Mathematica [77]. 
This gives 

q f[i] _ n /" (0!_1) sinh ( Lx p)[ x p( L (ci + c 2 ) + 1) - cic 2 ] + Lk p {c\C 2 + x|) cosh(Lx p ) 
2 J p 2x 2 [(cic 2 + x 2 ) sinh(Lx p ) + x p (ci + c 2 ) cosh(Lx p )] 

i-(d-i) 

/ [LI h (p,r) + I s (p,r,c 1 ) + I s (p,r,c 2 ) + I res (L;p,f,c 1 ,c 2 )} (A.24) 

J T> 



n 
2 



with the bulk, surface and L-dependent terms 

h(p,r) = ^-, (A.25) 
zx p 



and 



c\c 2 {c\c 2 L + ci + c 2 ) + Lxp - x 2 [(c 2 + c\)L + ci + c 2 ] 



Ires(L;p, T, Ci, c 2 ) — D „ „„ (A. 27) 

x p (ci + x p ) 2 (c 2 + x p ) 2 e 2L ^ - x p (x 2 - cf)(x 2 - c 2 .) 

of the integrand. Indefinite integrals Jb,s,rcs satisfying df Jb,s,res = -fb,s,res can be 
determined in a straightforward fashion. One finds 

J h (p, f) = Xp, J s (p, f, Cj) = In [x~ 1/2 (cj + Xp)], (A.28) 

and 



J res (L;p,f,Ci,c 2 ) = In 



(x p - ci)(xp - c 2 ) e" 
(ci + Xp)(c 2 + Xp) 



-2L> 



(A.29) 



Since reduces to p for f = 0, the integral J Jb would yield the usual bulk 
divergence oc A d if we regularised by cutting off the integral at p — A. It vanishes 
in the dimensional regularisation scheme because integrals of powers of p are zero. 
Likewise, f^ d ^ J s (p, 0,0) would yield a surface divergence cx A '" 1 (mod In A) if cut- 
off regularisation were used, but vanishes in our dimensionally regularised theory where 
only the Cj-dependent part contributes. This gives 

"(d-l) 

J s {p, 0, tj) = -tf d _i — — cj" 1 , (A.30) 



/ 



(d - 1) sin(Trd) 

where if<i_i is defined in equation (3.36). 

A straightforward combination of the foregoing results gives back equation (A. 22). 
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Appendix B. Computation of two-loop free-energy terms 



Starting from equations (3.29) and (3.31), we first consider the contribution arising from 
the part proportional to S mitm2 of A miimi)m2jm2 (Ci,C 2 ). This involves the integral 

' -^-i(^ + r)^ 2 , (B.l) 



p 2 + kl l + r 

where A d is defined in equation (3.64). From the result we see that a first series we 
must evaluate is J2m=i ^(^m + b 2 ) d ~ 3 . This can be done along lines analogous to those 
followed in Appendix A.l. Since the coefficient A m can be expressed as \ m = 2/gc lt c 2 ( K m) 
according to equations (2.19) and (A. 7), the function gci,c 2 drops out of the integrand 
of the contour integral. Thus the analogues of equations (A. 10) and (A. 17) become 



5> 

m=l 



2\d-3 



b 2 ) 



Res 

K=0 



1 

2vri 

b 2(d-3) 



dn 



73 



(^ + b 2 ) d - 3 N Cl , C2 (K) 



:«) 



i + cr 1 + c 2 - x 



+ 4 sin(7rrf) 



dt 



7T 



(t 2 



dn 



IT 



'2\d-3 



2\d-3 



+ b Z ) 

t + Ci t + C 2 
t-dt-C 2 



(B.2) 



.2* 



The integral over k exists when b > and d < 5/2. It can be computed. The result 

— (K 2 + b 2 ) d - 3 

7T 

4 



r(5/2-d) feM _ 5 



o 



(B.3) 



vrr(3 - d) 

provides the desired analytic continuation to values d > 5/2. For values of d = 4 — e 
close to the upper critical dimension d* = 4, both the term oc & 2 ( d_3 ) in equation (B.2) 
and the integral (B.3) approach zero as b — > 0. Hence 

-i -l 



2(d-3) 
m 



4sin(7r 



t + Ci t + C2 2t j 

t-d t-c 2 e 



(B.4) 



m=l 

for such values of d. 

A glance at equations (3.29) and (3.31) shows that the two-loop term of the free- 
energy density also involves the series 

* ^(^ + ^-3)/2 

with a = 0, 1,2. These can be evaluated in the same manner as the series (A. 4) and 
(B.2). We therefore simply quote the result one obtains after taking the limit b — > 0. It 

reads 

,.2<j+d-3 



5> 

m=l 



(^m + CD^l + C 2 ) 



2X, 



(d,a) 



2Y r 



(d,a) 



Ci,C 2 ' 



(B.5) 



where xj?'£) an< ^ ^c^A are the quantities defined by equations (3.38) and (3.39), 
respectively. 

From the above results the expressions for the bulk, surface and residual free-energy 
densities at f = given in equations (3.33), (3.35) and (3.42), respectively, follow in a 
straightforward fashion. 
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Appendix C. Calculation of the effective two-point vertex 7^) 

In this appendix we compute the eigenvalue k\ and the shift Stl(6i,c 2 ) for arbitrary 
nonnegative values of b\ and c 2 and derive the result (3.67) for the scaling function 
f2(ci,c 2 ). According to equations (3.57) and (3.58) the effective two-point vertex 7^ 
introduced in equation (3.65) is approximately given by 



k\ + St l = L~ 



2 n + 2u Tf 
+ — 2 L E 



P 2 + K 2 m 



(C.l) 



m=2 J P 

where the eigenvalues K 2 m are functions of the dimensionless bare surface enhancement 
variables C\ and C 2 . The momentum integration is straightforward. To compute the 
series Ylm=2 we substitute expression (3.31) for A 1)1;m)m , add and subtract the (m = 0)- 
term and then use the summation formula derived in Appendix A.l. We thus obtain 

n + 2u 



L 2 {kj + Sr L ) = «? 



„d-3 



3 2 
n + 2u 



= Ki 



3 8 



+ Ai«?" 3 -4A 1)lil)1 (Ci,C2)«? 



^4*-! I X] ^l,l,m,m(Ci,C2)/C^" 1 - Ai j i j i ) i(Ci,C 2 )/cf 
\m=l 



<2-3 



(C.2) 



The right-hand side has a logarithmic UV singularity at e = (simple pole) 
originating from the contributions proportional to the functions -X^iA an< ^ ^Ci A ■ Upon 
expressing -u in terms of the renormalized coupling constant u, its Laurent expansion 
becomes 



rhs(C2) 



n + 2u 



6 



E 



Ai«i 



2cr 



3(^.1) 



+ 0(e°) 



n + 2 (C 1 + C 2 )^(dC 2 + ^) 
^T MAl {C 2 + k 2 ){C 2 + k 2 ) + ° (e } - 



(C.3) 



Its pole term gets cancelled by the counter-term provided by the contribution oc u/e to 
k 2 (Ci,C 2 ) = kI(Lci, Lc 2 ) = k\{L[xZ c ci, LfiZ c c 2 ). To show this we write 



k\{LhZ c ci, Ljj,Z c c 2 ) = ref (L//C1, Ljj,c 2 ) + 2ki 



n + 2 



3e z — f r " J d(LfiCj) 
and calculate the partial derivatives using the implicit function theorem. This gives 



3=1 



dK-y{L[XCi, LflC 2 ) 



(C.4) 



E c > 

J'=l 



g«i(Ci,C 2 

0& 



E c - 



dR Cl ,c 2 (n)/dK 



K = K\ 



(Ci + C 2 )«i cos Ki + 2C\C 2 sin 



(2 + d + C 2 )«i sin ^ - (d + C 2 + C X C 2 - k\) cos k x 
(Ci+C 2 )ki(CiC 2 + «?) 



(C.5) 



CiC 2 (Ci + C 2 + C1C2) + (Ci + C 2 + C? + C 2 2 )/t? + «1 
Ai(Ci+C 2 )«i(CiC2 



2(C 1 2 + K f)(C 2 2 + ^ 
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Hence to the order of our calculation, k\ + Stl is UV-finite when expressed in terms 
of renormalized variables. From the implied UV-finite part of (C.2) the result for the 
scaling function fi(ci,c 2 ) given in equation (3.67) follows in a straightforward fashion. 

Appendix D. Behaviour of the shift 5r L (Ci,C 2 ) for small Cj 

In this appendix we calculate the shift 5fL(Ci,C 2 ) to linear order in C\ and C 2 . We 
start from equation (3.57), use equation (3.31) for A lilmm , substitute the small-Cj 
expressions (2.24) for k\ and ft^ >2 and expand in Cj. This gives 

Ai 5 i imim (Ci,C 2 ) _ I (Ci+C 2 ^ ' 1 



p 2 + ^(Ci, C 2 ) p 2 + (m — l) 2 7r 2 ^ 2(m — l) 2 7r 2 p 2 + (m — l) 2 7r 2 

2 1 ( D ' 1 ) 

+ o + o(Ci,c 2 ) L 

[p 2 + (m - l) 2 vr 2 ] 2 V V 

To determine the shift from these results, we interchange the integration over p in 
equation (3.31) with the summation over m. The required momentum integrals are of 
the form 

(d-i) 1 f « d ~ 3 for s = 1, 

Upon summing over m we encounter the series 

r(3/2-d/2)r(d/2-l/2) 



£»" = ^ Co - d) = rf f;;_ r ' rf /9 - ? «' = ( d. 3) 



and 

oo 

^(mvr)^ 5 = vr d - 5 C(5-rf). (D.4) 

m=l 

From the above results, equations (3.72) and (3.73) follow in a straightforward manner. 
Appendix E. Calculation of — f to two-loop order 

In this appendix we calculate the the difference — f of free-energy densities to two- 
loop order. Its one- loop contribution is given in equation (3.78). To determine its 
two-loop term, we must compute the two graphs appearing inside square brackets on 
the right-hand side of equation (3.79). The rightmost one of these (red) can be read 
off from equations (3.61) and (3.62) by setting 5i~l = 0. Upon adding the contribution 



from the remaining graph q in square brackets, we arrive at 

-L e A 2 d _ l <.K 2 ^ d 3 ^A lil5l)1 (Ci, C 2 ) — 2nf 3 ^ n d m 3 A lilim/m (Ci, C 2 ) \ 

^ m=l > 



un(n + 2) 



4! 
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un{n + 2) 
4! 



7" e 4 2 



K d-3 



.d-3 



|2A U , 1 , 1 (C 1 ,C 2 )-A 1 ] 



E 

<T,p=0 



2(7 



p 



(C 1 2 + / t ?)(C 2 2 + ^) 



X&P) _|_ y Wp) 

Cl,C2 Cl,C2 



(E.l) 



where the definition of C,- should be recalled from equation (2.14). The result contains 
pole terms; its Laurent expansion reads 

2 



rhs(E.l) = 



n + 2u ^ A^f 7 



48tt e 
n + 2 u 



(c\ + a 



+ 0(e°) 



= — n- 



Ai 



(Ci+C 2 )«?(CiC 2 



(E.2) 



24tt e 1 (C? + k?)(C| 



+ 0(6°). 



The pole terms are cancelled by those provided by the one-loop term /I 1 ' — / upon 
expressing the bare variables Cj in terms of their renormalized analogues Cj via 
equations (2.32) and (2.34). Upon adding the one-loop term one arrives at the two- 
loop result for /^r e5 ,R i'*" 1 given in (3.80) in a straightforward fashion. 
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